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ABSTRACT 

We investigate the multiple stellar populations in one of peculiar globular 
clusters (GCs) M22 using new ground-based wide-held Ca by and HST/WFC3 
photometry with equivalent passbands, confirming our previous result that M22 
has a distinctive red-giant branch (RGB) split mainly due to difference in metal 
abundances. We also make use of radial velocity measurements of the large num¬ 
ber of cluster membership stars by other. Our main results are followings, (i) The 
RGB and the subgiant branch (SGB) number ratios show that the calcium-weak 
( Ca-w ) group is the dominant population of the cluster. However, an irrecon¬ 
cilable difference can be seen in the rather simple two horizontal branch (HB) 
classification by other, (ii) Each group has its own CN-CH anti-correlation. How¬ 
ever, the alleged CN-CH positive correlation is likely illusive, (iii) The location of 
the RGB bump of the calcium-strong (Ca-s) group is significantly fainter, which 
may pose a challenge to the helium enhancement scenario in the Ca-s group. 

(iv) The positions of the center are similar, (v) The Ca-w group is slightly more 
centrally concentrated, while the the Ca-s is more elongated at larger radii, (vi) 

The mean radial velocities for both groups are similar, but the Ca-s group has 
a larger velocity dispersion, (vii) The Ca-s group rotates faster. The plausible 
scenario for the formation of M22 is that it formed via a merger of two GCs in a 
dwarf galaxy environment and accreted later to our Galaxy. 


1 Based on observations made with the Cerro Tololo Inter-American Observatory (CTIO) lrn and 0.9m 
telescopes, which are operated by the SMARTS consortium. 

2 Partially based on observations made with the NASA/ESA Hubble Space Telescope , obtained at the 
Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555. 
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Subject headings: globular clusters: individual (M22: NGC 6656) — Hertzsprung- 
Russell diagram - stars: abundances - stars: evolution 


INTRODUCTION 


The last decade has witnessed a drastic paradigm shift on the true nature of Galactic 
GCs. In particular, the conventional wisdom of the simple stellar population of GCs is no 


longer valid in almost all GCs in our Galaxy. In this regard, variations in 


in GC RGB stars have been known for more than three decades (1 Cohen! 19781) and the 


Na-0 anti-correlati on, for example, seen in almost all normal GCs f Garretta et al. 20091: 


igliter elements 


D’ercole et ah 2 008 ) is now generally explained by two star formation episodes (see also, 


Lee 


2010 1). in which the proto-stellar clouds of the second generation of stars must have been 
polluted by the gaseous ejecta from the first ge neration of the intermediate-mass asymptotic 
gia nt branch (AGB) stars (lD’Antona fc Ventura 2007) or the fast rotating massive stars 
( Decressi n e t al . 2007). 


Contrary to normal GCs showing th e lighter elemental abundance variations , pecu¬ 
liar GCs such as u Cen ( Lee et al. 1999 : Bedin et ah 2004 ) and M22 ( Lee et al. 2009a : 
Marino et al.l 2009. 2011) show spreads in the heavy elemental abundances which synthesized 
and supplied by the core-collapsed supernovae. To retain ejecta from energetic supern ova 
explo sions requires much more massive systems than typical Galactic GCs (iB aumgar dt et al. 
20081 ). Furthermore, some_of massmyGCs in our Galaxy appear to have different kinematics 
than the normal GCs ( Lee et all 2007). As a consequence, these peculiar GCs are thought 
to h ave an extrag alactic origin, such as being the remnant of dwarf galaxies (for example, 


see 


Bckki & Freeman! 2 0031 ). and they are building blocks of our Galaxy. 


One of the early hints on the peculiarity of M22 can be found in iHesser et al.l (119771) 


and Hess e r & Harris (1979). They studied M22 using the DDO photometry and the low 


resolution spectro scopy, finding heterogeneous light elemental abundances in M22 RGB stars. 
Hesser & Harris also n oted t hat M22 may share the similar elemental abundance anomalies 


of c o Cen. Later Norr is fe Freeman (119831) found a variation in calcium abundance as well 


as those of the CN band and the G band strength in M22 RGB stars. They suggested that 
calcium abundances among M22 RGB stars can vary up to A[Ca/H] « 0.3 dex and it is 
positively correlated with the carbon abundance. 

Regarding variations of heavy elemental abundances including calcium in M22, more 


intriguing results have been emerged in recent years. In our previous work, Lee et ah (2009a) 
studied M22 using the extended Stromgren photometry and we found the distinctive RGB 
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split in the hk index [= (Ca — b) — (b — y)}. The hk in dex o f RGB stars is a measure 


of the Ca II H & K absorption strength (jAnthonv-Twarog et al.l Il99ll) and ou r previous 


result indicated discrete calcium abundances among M22 RGB stars (see also Lim et ah 
20151 ). We also showed that two discrete RGB sequences have indeed different iron and 


calcium abunda nces by using elemental abundance measurement fr om the high resolution 
spectroscopy bv lBrown fe Wallersteinl (119921 1 . iDa Costa et al.l (120091) also obtained the same 
result from observations of the infrare d Ca II triplet that M22 RGB stars show a substantial 
intrinsic metallicity spread. Finally, Marino et ah (l2009l 2011) performed high resolution 
spectroscopic studies of RGB stars in M22. In addition to the _spread in metallicity and the 
calcium abundance among RGB stars in M22, Marino et ah (boil ) found that the metal-rich 
group of stars have enhanced CNO, which may account for the the split in the SGB, and 
s- process-elemerit abundances. Furthermore, each stellar group in M22 appears to exhibit 
its own Na-0 anti-correlation. Spectroscopic studies of the M22 HB stars also showed that 
the the true nature o f HB populations m ay be as complex as that of RGB populations (e.g. 
Marino et al.l i2014at Grafton e t ah 2014). All the available observational evidences up to 
date may suggest that M22 must have retained not only the material from SNe but also 
from low-mass and intermediate-mass AGB^, reflecting that the formation history of M22 
must have been very complex compared to those of the bulk of normal GCs in our Galaxy. 

This is the first of the series of papers addressing the multiple stellar populations of 
Galactic GCs based on the homogenous Ca by photometry. In this study we present a new 
wide-fiel d Ca by photometry of M22. As we have extensively demonstrated in our previous 
work (jLee et al.l 2009af) . the narrow-band photometric study of GCs can provide an efficient 
means to investigate multiple stellar populations in the Galactic GC systems. Here we 
present more evidences that make M22 one of peculiar GCs in our Galaxy. 


2. SEJONG Ca by SURVEY 

Sejong Ca by survey program has been launched in July, 2006 after two years of prepa¬ 
ration. As an official major partner of the SMARTS consortium, we acquired the guaranteed 
access to the small telescopes operated at the CTIO since 2006. For 8 years, we used about 
270 nights in total, about 120 of which were photometric, for this survey program. 

The main goals of our survey program are three fold. Using the extended Stromgren 
system (Stromgren by filters plus the Ca filter), we aimed at obtaining followings; i) Homo- 


1 See also Decressin et ajJ (2007) for the fast rotating massive stars (FRMSs) as a source of the Na-0 


anti-correlation of GCs. 


































M22 


4 


geneous photometry of individual GCs without prerequisite information of the cluster. We 
tried to maintain a flux limited blind survey for Galactic GCs. The main purpose of our GC 
survey program is to study multiple stellar populations of Galactic GC systems to shed more 
light on the formation of the GC systems and furthermore on the formation of our Galaxy, 
ii) Metallicity distribution of selected fields in the Galactic bulge, in particular we are very 
interested in the metal-poor regime of the Galactic bulge, iii) Metallicity distributions and 
multiple stellar populations of nearby dwarf galaxies. 


Our survey program is con sisted of three phases. The Phase I was from July, 2006 to 
November, 2010 (see Lee et al.l 2009a,b). During this period, we used Stromgren uvby and 
Ca Liters provided by the CTIO. The CN band at A 3885 A is often very strong in the RGB 
stars and the CTIO Ca filter was originally designed to avoid the CN band contamination. 
However, it was suspected that the CTIO Ca Liter had undergone degradation due to aging 
effect and the original transmission function had been altered to the shorter wavelength. As 
a result, the CTIO Ca Liter is suspected to be suffered from the CN band contamination. 

Since March 2011 (Phase II), we used our own Stromgren by and new Ca Liters, all of 
which were manufactured by Asahi Spectra, Japan. Our new Ca Liter was carefully designed 
to have very similar Liter band width and pivot wavelength as those of F395N Liter in the 
Wide-Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST). In Figure [U we 


show a comparison of Liter transmission functions between that in Anthony-Twarog et al. 


(119911 ) and that of our new Ca Liter measured with collimated beam by the manufacturer of 
the Liter. Both Liter have similar full-width at the half maximums (FWHMs), approximately 
90 A but our new Ca Liter has a more uniform and high transmission across the passband, 
dropping more rapidly at both edg es. As s hown in the Lgure^ the CN band at A 3885 A lies 
on the lower tail of the Ca Liter by Anthon y- Twaro g et al . (11991 1 while our new Ca Liter is 
designed to be completely free from the CN band contamination. 


In Figure [H we also show the Liter transmission functions for the HST WFC3 photome¬ 
try system. The HST WFC3 F547M Liter is supposed to be equivalent to the ground-based 
Stromgren y Liter. However, the bandwidth of the F547M is twice as broad as the conven¬ 
tional Stromgren y Liter. 


The Phase III started in April 2013, we added a new Liter system, JWL39 , which allows 
us to measure the CN band at A 3885 A in combination with our new Ca Liter. Note tha t 
the bandwidth of JWL39 is similar to that of DD038 ( McClure fc van den Bergh 1968 1. 
but the pivot wavelength of the Liter is at ~ 3900 A. We show the transmission function of 
our JWL39 Liter in Figure |U 


In our current work, we will present Ca by photometry of M22 taken during the Phases 
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II and III only. Wc will present the results using the new JWL39 filter system in forthcoming 
papers. 

There are at least two elements of strength in our survey program. First, we used 
the same instrumental setups, the CTIO lm telescope and the Y4KCam, for most of our 
observations except for two runs in 2012, due to the maintenance of the Y4KCam. Second, 
almost all the data have been obtained by the author of the paper and all the data have 
been reduced and analyzed by the author of the paper, which guarantees homogeneity of the 
final results. 


3. OBSERVATIONS AND DATA REDUCTION 
3.1. Observations 

New Ca by observations for M22 were made in 44 nights, 23 of which were photometric, 
in 10 runs from March 2011 to May 2014 using the CTIO 1.0m and 0.9m telescopes. The 
CTIO 1.0m telescope was equipped with a STA 4k x 4k CCD camera, providing a plate 
scale of 0.289 arcsec pixeD 1 and a field of view of about 20 x 20 arcmin. The CTIO 0.9m 
telescope was equipped with the Tektronix 2048 No. 3 CCD, providing a plate scale of 0.40 
arcsec pixel” 1 and a field of view of 13.5 x 13.5 arcmin. For most of our observations, we 
used the CTIO 1.0m telescope except for one semester. The Y4KCam mounted on the CTIO 
1.0m telescope had been under inspection at Steward Observatory, Arizona, USA during the 
spring semester of 2012, and we used the CTIO 0.9m telescope for two runs in April and 
July 2012. 

The total exposure times for M22 were 90,900 s, 20,930 s and 9,145 s for our new Ca, 
Stromgren b and y , respectively. We show the journal of observations in Table [0 


3.2. Data Reduction 

The raw data were processed through the standard IRAF 2 packages using twilight or 
dawn sky flat images to calibrate science exposures. Since both CCD cameras are equipped 
with an iris type shutter, the illumination across the CCD chips is not uniform and the 


2 IRAF (Image Reduction and Analysis Facility) is distributed by the NOAO, which are operated by 
the Association of Universities for Research in Astronomy, Inc., under contract with the National Science 
Foundation. 
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shutter shading correction is applied for the p hoto metric standard frames with the exposure 
times shorter than 10 s (see 


Lee fe Pogge 20151 . for detailed discussion). Note that the 
difference in the integration time between the center and the edge o f th e CC D is 7C 


for t he Y4KCam. The most of photome tric standards bv lTwarog fe Anthonv-Twarog 


msec 


1995) 


and lAnthonv-Twarog fe Twarogl (119981) are too bright even for the small telescope like the 
CTIO 1.0m telescope. For most cases, we intended to have exposure times longer than 10 s. 
Although defocusing observing technique has been frequently used to observe isolated bright 
stars in the held, such as photometric standards, the most of the standard stars are still too 
bright to be observed with a integration time of longer than !=s 10 s. Without any viable 
options, science frames with a short integration time are often necessary in order to secure 
wide range in colors of photometric standards, especially in y passband. 

The held of view of the hnal combined science image of M22 is about 55 x 55 arcmin, 


which is more than 7 times larger than that of our previous study of the cluster dLee et al. 
2009aj). The photometry of the cluster and photometric standard frames were analyzed 
using DAOPHOTII, D AOGROW, ALLSTAR and ALLFRAME, and COLLECT-CCDAVE- 
NEWTRIA L packages (Stetson 1987 , 


1998, 


described in Lee & C ar ney (1199911 and Lee et ah (2014). In order to derive the transformation 


1994 . Il995l; iTurnerl 1199511 following the method 


relations from the instrumental system to the standard system, we adopt the following 
equations: 


Uinst = VSTD + OCy(b — y)sTD + fiyX + 'fy, ( 1 ) 

binst = bsTD + &b(b — y)sTD + AW + 

Gainst = CasTD + OlCaib — y)sTD + PcaX + ShksTD + 7 Ca, 


where X is the airmass and ysTD, bsTD, CasrD, (b — y)sTD and hk STD are the magnitudes 
and color on the standard system. The total number of stars measured in our M22 held from 
our ALLFRAME run was more than 518,000. Most of stars measured in our held turned 
out to be Galactic bulge stars. 


In Figure E] and Table [21 we show compa rison s of our photpmetry with those of pre¬ 
vious studies bv lAnthonv-Twarog et al.l (1199511 and iRichter et al.l (1199911 . In the hgure, the 
differences are given in the sense of other works minus our work. It can be seen that our 
results are in good agreement with previous studies within measurement errors. It should be 
emphasized th at total exposure times o f our obser vations in each passb and are much greater 
than those of lAnthonv-Twarog et al.l (119951) and IRichter et al.l (1199911 . and, therefore, the 
photometric quality of our work is believed to be superior to those of previous works. 


Finally, astrometric solutions for individual stars in our held have been derived using 
the data extracted from the Naval Observatory Merged Astrometric Dataset (NOMAD, 
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Zachairias et ah 2004j). We achieved the rms scatter in the residuals of less than 0.07 arcsec 


using the IRAF IMCOORS package. Then the astrometric solution was applied to calculate 
the equatorial coordinates for all stars measured in our science frames. 


3.3. Completeness tests 


To examine the completeness of our photometry as a function of V m agnit u de, w e 
performed a series of artificial star experiments (see for example, Stetso n & H a rris 1988)). 
We selected two square regions with a size of about 15 x 15 arcmin from the central part and 
the outer part of the cluster. We constructed a fortran program to distribute 1,000 stars in 
both fields by adopting the observed radial stellar number density and the luminosity profiles 
in our study. Note that the number of artificial stars used in our experiments corresponds 
to about 1.3% and 5.6% of the total measured number of stars in the central and the outer 
part of the cluster, respectively. This number of artificial stars added to the observed images 
is carefully chosen not to dramatically change the crowdedness characteristics between our 
data reduction procedures and our artificial star experiments. 


We used the observed point-spread functions (PSFs) and the DAOPHOT’s ADDSTAR 
task to generate 20 artificial star images in each field. Following the same data reduction 
procedure that we described above, we derived the completeness fractions as a function of 
V magnitude. We show our results in Figure [3j As shown in the figure our experiments 
suggest that our photometry is complete down to V ~ 18.5 mag both in the central and in 
the outer field of the cluster for the ground-based observations. As will be discussed below, 
it should be emphasized that incompleteness at the fainter magnitude regime does not affect 
our results presented in this work. 


4. RESULTS 

4.1. Color-Magnitude Diagrams 

4-1.1. Ground-based, observations: field star contamination 

In Tabled we provide our photometric data for stars brighter than V — 19 mag. Figured 
shows the color-magnitude diagrams (CMDs) of bright stars in the M22 field. Note that, 
for the sake of clarity, we show stars within 10 arcmin from the center of the cluster. As 
can be seen in the figure, the double RGB sequences of M22 can be clearly seen, especially 
the distinctive split in the RGB stars in the hk versus V CMD. Note that the ml versus V 








M22 


CMD is from our previous study of the cluster ([Lee et ah 2009a) taken during the Phase I 
period. 

Since M22 is located towards the Galactic bulge (l = 9.9°, b = —7.6°), the contribution 


from the field star contamination is naturally expected to be very large (see iLee et ah 2009a. 


for detailed discussions). We attempted to remove the off-cluster field stars in our bright RGB 
sample by using multi-color CMDs following the method described in ILee et aL (2009a. b), 
which turned out to be very successful as discussed below. 

In Figure [5]-(a) and (b), we show CMDs of the M22 field using all the stars measured 
from V = 11.5 mag to 17.5 mag. In the figures, blue solid lines indicate the RGB boundary 
grids in the (b — y ) versus V and in the hk versus V CMDs. In order to prepare these RGB 
boundary grids, we made use of the fiducial sequence of the model isochrones by Joo & Lee 
(120131 1. which were kindly provided by S.-J. Joo, an d we a dopted the distance modulus and 
the foreground reddening value for M22 by Harris (119961 ). Then we added small amount 
of color offset values, A(b — y) ~ ± 0.04 mag and A hk ~ ± 0.07 mag, in the transformed 
model isochrone sequences to make the final RGB boundary grids. Although our definitions 
of the RGB boundary grids for (b — y) and hk indexes appear to be arbitrary, our approach 
appears to work with great satisfaction as shown below. 

The total number of RGB stars detected in the enclosed region in Figure [5]- (c), is 3571. 
When plotted in the hk versus V plane, it can be clearly seen that about 42% of RGB stars 
selected in the (b — y) versus V plane does not belong to M22, i.e. either the disk or the bulge 
po pulations wi th large hk index values. A comparison with the radial velocity measurements 
by lLane et al.l (120091 ) shows that 12 stars out of 307 stars in common turn out to be radial 
velocity non-member stars. If we take this as a face value, the fraction of the inclusion of 
the off-cluster field population in our M22 RGB selection procedure is 3.9 ± 1.4%, which 
appears to be too low to dramatically change the results presented in this paper. Without 
12 radial velocity non-member stars, we have, 2,054 RGB stars, from V = 12.0 mag to 16.5 
mag, in total. 


The validity of our RGB selection procedure can be confirmed by the surface-brightness 
profile (SBP) of the RGB stars in M22. In Figure El- (d), we show the SBP of the red-clump 
stars in the Galactic bulge [the stars lie inside a cyan box in Figure [5]- (a)]. As expected, 
the SBP of the bulge population is almost flat against the radial distance from the center of 
M22. Note that the mild fluctuation seen in the small r region, log(r/arcsec) < 2, is due to 
the small number of the bulge red-clump stars (n < 15) detected in this region. 

Figure [6]- (a)-(c) show the SBPs of the Ca-w group, the Ca-s group and all RGB stars in 
M22 (see below for details on the definition of the Ca-w and the Ca-s groups). With slightly 
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different zero-point offset values for each group of stars, our SBPs for M22 RGB stars are in 
excellent agreement with the Chebyshev polynomial fit of M22 by Trager et ah (1199511 up to 
more than 10 3 arcsec from the center of the cluster. Our investigation of the SBPs of M22 
RGB stars strongly suggests that the contamination of the off-cluster held star in our result 
is not severe, consistent with our previous result from the comparison of the radial velocity 
measurements of M22 RGB stars. 


4-1.2. Ground-based observations: Ca-w versus Ca-s 


As shown in Figures [I] and 0 M22 exhibits the distinctive double RGB sequences in the 
hk versus V CMD, which is mainly due to the differ e nce in the heavy elementa l abundances 
between the tw o groups (see als o iLee et al.l l2009al; iMarino et al.l 120091 . 1201 111 . As in our 
previous works (ILee et ahl 2009al lbl). we define the calcium-weak (Ca-w) group of stars with 
smaller hk index and the calcium-strong (Ca-s) group of stars with larger hk index at a 
given V magnitude. 


In the left panel of Figured we show 2,054 RGB stars from Figure [5]-(d) on the A hk 
versus V plane, where A hk is defined to be the difference in the hk index of the individual 
RGB stars with respect to the hk index of the fiducial sequence of the Ca-w RGB group. 
As shown, the differences in the hk index between the mean values of the Ca-w and the 
Ca-s depend slightly on V magnitude, which is naturally expected if the difference in the 
absorption strengths of Ca II H and K lines is mainly responsible for the separation in the hk 
index of the RGB stars. Due to this magnitude (more precisely, surface gravity) dependency 
on the A hk, we divide M22 RGB stars into four different magnitude bins and we show 
histograms of the A hk distributions of RGB stars of each magnitude bin in the right panels 
of Figured all of which appear to show double peaks. In order to separate the two different 
groups of RGB stars based on the A hk distribution, we adopt the expectation maximization 
(EM) algorithm for the two-component Gaussian mixture distribution model. In an iterative 
manner, we derive the probability of individual RGB stars for being the Ca-w and the Ca-s 
groups. In the left panel of Figured stars with P(w\xi) > 0.5 from the EM estimator are 
denoted with blue dots which corresponds to the Ca-w group, while P(s\xi) > 0.5 with red 
dots which corresponds to the Ca-s group. In the right panel of the figure, we show the 
distributions of the Ca-w and the Ca-s groups with blue and red solid lines, respectively, 
and we show the number ratio n(Ca-w)m(Ca-s) for each bin. In total the number ratio 
between the Ca-w and the Ca-s groups is n(Ca-w):n(Ca-s) = 68.8:31.2 (±3.3). 


Our observed number ratio between the two groups does not necessarily reflect the initial 
stellar number ratio or the initial mass ratio between the two groups, since both groups 
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have different initial chemical compo sitions, different masses at a given luminosity, and 


furth ermore slightly different ages (e.g. JLee et al.ll2009atiMarino et al.ll2009l.l201ll: I Joo fe Lee 
2013) . In order to take care of the heterogeneous evolutionary effects between the two 


groups, we construc ted a n evolutionary population synthesis model using the theoretical 
model isochrones by Joo fe Lee ( 20131 ). We populated 10 7 artificial stars in each group using 
the Salpeter’s initial mass function (IMF) and we generated 50 different sets for each group. 
We, then, compared the number ratios between the two groups in the magnitude range of 
12.0 < V < 16.5 mag, finding n(Ca-w)\n(Ca-s) = 48.3:51.7 (± 0.3), in the sense that RGB 
stars in the Ca-s group are slightly more numerous than those in the Ca-w group at a fixed 
total number of stars in both groups. Note that using other t ypes of IMF does not affect 
our results presented here, which was also pointed out by [Joo_& Lee (2013). As an exercise, 
we calculated the number ratio using the universal IMF by Kroupa (120011 ) and we obtained 
n(Ca-w):n(Ca-s ) = 48.4:51.6 (± 0.3), in excellent agreement with that from the Salpeter’s 
IMF within the 0.1% level. 

Figure [8] shows the mean theoretical luminosity functions (LFs) for the Ca-w and the Ca- 
s groups. Note that the abscissa of the figure is M v and the grey shaded area is corresponding 
to M22 RGB stars with 12.0 < V < 16.5 mag. As can be seen, the theoretical LFs for the RGB 
stars in both groups are very similar. If we consider the evolutionary effect, the intrinsic or 
the true number ratio between the two groups becomes n(Ca-w):n(Ca-s ) = 70.2:29.8 (±3.3). 
Note that our new RGB number ratio is in good agreement with our previous rough estimate 
of the number ratio between the two groups, ~ 0.70:0.30 ( Lee et al. 2009al) . 


As discussed by iLee et al.l (]2009al ) and IMarino et al.l (120091. 120111 ). the double RGB 
sequences in M22 is due to the difference in chemical abundances between the Ca-w and 
the Ca-s groups, in particular the difference in the calcium abundance measured with the 
Ca filter and in heavy e lements, which was fir s t hint ed by the high-resolution spectroscopy 


of sev e n RG B stars by iBrown fe Wallersteinl (119921) and confirmed later by IMarino et al. 
(120091 . 2 0111) with the expanded sample of RGB stars of the cluster. More detailed dis¬ 


cussion on the difference in the chemical abundances between the two groups will be given 
below. Also worth emphasizing is the fact that the differential reddening effect and the 
contamination from thepoff-cluster populations cannot explain the double RGB sequences 
in M22 (iLee et al.l l2009al) . Assuming that the c onventional interstellar r edden ing law for 
the extended Stromgren photometry system by A nthonv - Twarog et ah (J199ll ), E(hk ) = 
—0.16 xE(b — y) ~ — 0.12x77(5 — V), the mean hk difference between the two RGB groups, 
A hk = 0.09 mag, can be translated into A (b — y) — —0.56 mag, A (B — V) = —0.75 mag 
and AV = —2.33 mag, in the sense that the metal-rich Ca-s RGB stars becomes bluer and 
brighter due to the foreground reddening. Indeed, this is not the case. 
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4-1.3. Ground-based observations: ml versus V 


The Stromgren v passband includes the CN band at A 4215 A and, as a consequence, 
ml index [= (v — b) — (b — y)\ is capable of distinguish the difference in the CN absorption 
strengths. In fact, the bimodal distribution in the ml index of M22 RGB stars was known for 
decades and it has been frequently attributed to the bimodal distribution in the CN abun¬ 


dances (iNorris fe Freeman Ill983l: lAnthonv-Twarog et al.lll995t iRichter et al.lll999t iLee et al. 


2009aJ). It is widely accepted that the variation of lighter elemental abundances, such as C, 
N and O, is thought to be res ulted from chemical pollution by the intermediate-mass AGB 
stars (see, for example, Ventur a et al. 20011). The potential trouble with the ml index is that 
the split between the double populations at the lower RGB sequence is not readily distin¬ 
guishable in the ml versus V CMD, most likely due to the strong surface gravity sensitivity 
on the CN band strength, in the sense that th e CN band strength is i nversely correlated w ith 
the surface gravity of stars (see for example, iTripicco fe Belli Il99ll: ICannon et al .1119981) . as 
shown in Figure^ 


On the other hand, as Lee et ah (j20 09af ) discussed extensively, the hk index is a measure 
of calcium abundances of individual stars and, therefore, the difference in the hk index at a 
given luminosity indicates the heterogeneous calcium and other heavy elemental abundances 
between the double RGB populations in M22 , which is confirmed by the high-resolution spec¬ 


trosc opic study of RGB stars in the cluster flBrown fc Wallersteinlll992l; iMarino et al.l 12009 


20111). Calcium and heavy elements can only be supplied through supernovae explosions of 
massive stars and the origin of the split in the hk index is intrinsically different from that in 
the ml index in M22, if the ml index is the most sensitively depends on the CN abundances 
as previous investigators perceived frequently. It is important to note that, however, since 
the Stromgren v passband includes numerous metallic lines, the ml index also sensitively 
behaves in the difference in the heavy metal abundances. 

In the top panels of Figure [9l we show the CMDs of M22 again. For the ml versus 
V CMD, we show stars with (b — y) > 0.45 mag to prevent the degeneracy of the location 
of the HB and the RGB stars in the ml versus V CMD. We also show the Padova model 
isochrones with [Fe/H] = —1.8 dex, comparable to that of the Ca-w group, and —1.5 dex, 
comparable to or slightly more metal rich than that of the Ca-s group, with a fixed [CNO/Fe] 
ratio OBressan et al, 2012). From the figure, it is palpable that the bifurcation in the RGB 
sequences on the ml versus V CMD is mainly due to the difference in the metal abundance. 
In addition, the higher mean [CNO/Fe] ratio in the Ca-s groupie may intensify the difference 
in the ml index between the two groups. 


2 The s-process-rich group 
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A comparison of M22 with two Galactic GCs, M55 and NGC 6752 (Lee 2015 in prepa¬ 
ration), may help to elucidate the nature of the double RGB sequences. Note that the 
metallicity of M55 is [Fe/H] = —1.9 dex, equivalent to the Ca-w group in M22, while that 
of NGC 6752 is —1.5 dex, equivalent to the Ca-s group in M22. It is well known that the 
HB morphologies of both clusters are different. M55 contains only the blue HB (BHB) pop¬ 
ulation, while NGC 6752 contains the BHB and the extreme blue HB (EBHB) populations. 
The CN distributi ons o f both cl uster s are also different. M55 appears to show an unimodal 


CN distrib ution (IKavser et a 


NGC 6752 (Norris et al. 1981). 


20081), in sharp contrast to the bimodal CN distribution of 


In the bottom panels of Figure[9l we show the composite CMDs of M55 and NGC 6752. 
In order to make these composite CMDs, we tweaked the color indexes and the visual mag¬ 
nitude of individual stars in NGC 6752 to match with those in M55; A (b — y) = 0.034 mag, 
Aral = —0.008 mag, A hk = —0.001 mag and AH = 0.654 mag. In the ml versus V CMD, 
we do not plot the HB stars in M55 and NGC 6752 to avoid the confusion with the RGB 
stars. We emphasize that these composite CMDs can successfully reproduce the photometric 
characteristics of M22. Especially, the CMDs of M22 and the composite GCs show remark¬ 
able resemblance in the HB morphology and the double RGB sequences both in the ml and 
the hk indexes, suggesting that the double RGB sequence in the ml index of M22 is likely 
due to the difference in the mean heavy metal abundance. If so, the origin of the RGB split, 
in the ml index is the same as that of the hk index in M22. 


Nonetheless, using the hk index as a metallicity indicator has advantages. As Anthony-Twarog et al. 


1 19911) noted, the hk index is about three times more sensitive to metallicity than the ml 
index is, for stars more metal-poor than the Sun, and it has half the sensitivity of the ml 
index to interstellar reddening. The annoying trouble with the ml index is that it loses 
sensitivity to metallicity changes in the regime below [Fe/H] « —2 dex, especially for stars 
with higher surface gravity. 


4-1.4- Hubble Space Telescope Observations 

We also carried out the HST WFC3 photometry for the central part of M22. In Fig¬ 
ure C01 we show the preliminary results from HST WFC3 observations of the cluster (PI : 
J.-W. Lee and PID : 12193) using the similar passbands that we used in our ground-based 
observations. 

The left panel of the figure shows the (m.F 467 M — ^f 547 m) versus ra.F 547 M CMD, which is 
equivalent to the ground-based (b — y ) versus V CMD, near the the subgiant branch (SGB) 
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and the faint RGB regions of M22. Similar to the gronnd-based (b — y ) versus V CMD, the 
split in the RGB or SGB can not be readily seen in the figure. 

On the other hand, the hk s tmag [ = ("IF395N - rn F467M ) - (m F467M - m F547M )] versus 
m F 5 47 M CMD shows the distinctive split in the the SGB and the RGB sequences of M22, con¬ 
firming the results from our new ground-based observations discussed above. Note that the 


SGB split in M22 has been already known in other HST passbands by Marino et ah (2009 


2012 ) and the prese nce of the SG B split is likely due to variations in the CNO abundance 
(iCassisi et al.l 120081) . Especially, iMarino et al.l (1201211 performed a low-resolution spectro¬ 
scopic study for a large number of SGB stars in M22, finding that the elemental abundance 
trend of the faint SBG (fSGB) group is consistent with that of the s-process-rich group in 
their early study (IMarino et ahl 2 01ll. i.e. Ca-s in this study), while the elemental abundance 
trend of the bright SGB (bSGB) is consistent with that of the s-process-poor group (i.e. the 
Ca-w group), demonstrating that the observed difference in the CNO abundance can account 
for the observed SGB split. 


Mil on e et al. (120091 1 devised a rather clever but complicated method to estimate the 


number ratio between the fSGB and the bSGB in NGC 1851. However, there are at least 
thre e aspe cts that make difficult for M22 to be analyzed using the method described by 
Milone et al.l (120091 ). First, the spatial stellar number density in a fixed area on the sky (i.e. 
the number of stars/d") of M22 is about 20 times smaller than that of NGC 1851. Therefore, 
the number of the SGB stars detected in our HST observations for M22 are much smaller 


than that in NGC 1851 by Milone et ah| (120091 ). The small number of stars makes it difficult 
to delineate a clear SGB fiducial sequence for M22. Secondly, the foreground interstellar 
reddening value of M22, E{B — V ) = 0.34 , is much larger than that of NGC 1851, E(B — V ) 
= 0.02 (lHarrislll996 ). Furthermore, there may exist the differential foreground reddening 
effect across M22, although the differential reddening is not the main reason for showing the 
double RGB or SGB sequences in M22 as we have already discussed. Lastly, the curvature 
of the SGB stars and the location of the transition region between the SGB and the base of 
the RGB are different between the HST hksTMAG and the F606W — F814IF photometric 
systems. The essential idea of the method proposed by Milone et al. (2009) is to estimate 
the intrinsic or the true number ratio between the two SGB populations, derived from the 
apparent number ratio by taking care of the evolutionary effect as we have done for the RGB 
stars. It also should be important to point out that their method is arbitrary and model 
dependent, as they have noted. 


Before turning to the observed number ratio between the two SGB groups in M22, we 
investigate the methodological aspect by employing the evolutionary population synthesis 
models that we constructed for the RGB stars above. In the left panels of Figure [HI we 
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show artificial CMDs around the SGB region of M22 for different photometric systems. In 
orde r to cons truct the evolutionary population synthesis models, we use the model isochrones 
by Joo & Lee (2013) with the Salpeter’s IMF. The photometric errors are estimated from 
the measurement errors around the SGB region in our HST observations. We populate 10 8 
artificial stars in each group. In the figure, blue dots denote the bSGB group while red dots 
the fSGB in M22. For the sake of clarity, we only show a small fraction of stars in the figure. 
We then count the numbers of SGB stars for each group at a fixed color range (the vertically 
shaded areas) and at a fixed magnitude range (the horizontally shaded areas) returned from 
our calculations. 


The middle panels of Figure fill show the distributions of SGB stars in the fixed color 
ranges for the different photometric systems, while the right panels show those in the fixed 
magnitude ranges. Similar to our RGB number ratio estimates, we calculate the differences 
in the magnitudes or in the color indexes against those of the fiducial sequences of the bSGB 
group. The first number ratios in each plot are those returned from the EM estimator 
assuming the two-component gaussian mixture model while the number ratios inside the 
parenthesis are the input values used in our evolutionary population models. Note that 
these number ratios in the parenthesis reflect the evolutionary effect between the two SGB 
groups at a given color or magnitude ranges. In the figure, our adopted ranges in each color 
index may be rather arbitrary but we intend to maximize the available number of stars 
maintaining the clear separation between the two SGB groups. 

First, the SGB number ratio returned from the EM estimator between the two groups 
at the fixed (F606W — F814W) color is n(bSGB):n(fSGB) = 42:58, while that of the input 
value from the evolutionary population synthesis models is 39:61. On the other hand, the 
SGB number ratio determined at the fixed F606W magnitude appears to be more uncertain. 
We obtained n(bSGB):n(fSGB) = 63:37 with the input value of 56:43. 

Next, the SGB number ratios determined from the b — y versus V CMD appear to be 
the least reliable in our simulations. We obtain n(bSGB):n(fSGB) = 44:56 for the fixed 
color and 72:28 for the fixed magnitude ranges, which number ratios deviate largely from 
the input values of 38:62 and 65:35, respectively. 

Lastly, the SGB number ratios determined from the hk versus V CMD appear to be the 
most reliable. Due to the combined effect of the mean heavy metal and the CNO abundances 
on the location of the SGB sequence on the hk versus V CMD, the separation in the hk index 
between the two SGB groups is the most conspicuous. The EM estimator securely retrieve 
the input values within an uncertainty of < ± 0.2. Since our simulations suggest that the 
total number of stars available in our simulations for the hk versus V CMD is about 3 times 
larger for the case with the fixed magnitude range, we estimate the observed SGB number 
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ratio from our HST observations of the cluster using the method with a fixed magnitude 
range, i.e. using the distribution of individual stars against A hk, as for the RGB stars. Also 
importantly, using the method with a fixed magnitude range has an advantage in correcting 
the evolutionary effect. The SGB number ratio of n(bSGB):n(fSGB) = 58:42 is, in fact, the 
evolutionary effect correction factor shown in Figured] (the yellow shaded area) and applying 
the correction factor of the differential evolutionary effect based on the V magnitude is more 
direct. 


In Figure H2]-(a), we show the /i/cstmag versus m F547 M CMD around the SGB region 
of M22. Also shown is the fiducial sequence for the bSGB group. Following the similar 
procedure that we used for the RGB stars, we obtain the number ratio of n(bSGB):n(fSGB) 
= 69:31 (± 6). If we apply the correction factor for the differential evolutionary effect 
between the two group, the intrinsic number density becomes n(bSGB):n(fSGB) = 62:38 (± 
6). This SGB number ratio is in agreement with the RGB number ratio, n(Ca-w):n(Ca-s ) 
= 70:30 (± 3.3), within the measurement errors. 


As mentioned above, Marino et al. (120091. see references therein) already discovered the 
SGB split in their HST Advanced Camera for Surveys Wide Field Channel (ACS/WFC) 
observations of M22 using F606W and F814W passbands, although the separation between 
the fSGB and the bSGB populations in the m F60 6w — wt-fsmw versus m F si 4 w CMD is not 
as clear as that in hfcsTMAG versus m F 547 M CMD as shown here. Their SGB number ratio is 
n(bSGB):n(fSGB) = 62:38 (± 5), which is in excellent agreement with our true SGB number 
ratio n(bSGB):n(fSGB) = 62:38 (± 6). 


Finally, it should be worth mentioning that the hk separation between the two SGB 
groups in the observed CMD (Figure flOl) is more ambiguous than that of the simulated 
CMD with a proper treatment of the photometric measurement errors (Figure fill). This 
may indicate that not only the differential reddening but also the metallicity spreads in 
the double SGB populations are necessary to perform more realistic simulations. In fact, 


Marino et al.l (2011]) reported rather large metallicity spreads for the two groups of RGB 
stars; 0.10 dex for the s-process-poor group and 0.05 dex for the s-process-rich group. As 
shown in Figure [HI the (b — y) width of M22 is compatible with that of the composite GC 
(i.e. NGC 6752 + M55), where the (b — y) locus of the RGB sequence is less sensitive to the 
variation in the heavy metal abundance, while the hk width of M22 is broader than that of 
the composite GC. Given the current interstellar reddening law (see below for more detail), 
the reasonable explanation for the broader RGB sequence of M22 in the hk index is likely 
due to the spread in the heavy metal abundances of both stellar groups in M22. 
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4.2. Differences in Chemical Abundances between Two RGB Groups 


Lee et ah (I2009al ) showed that the distinctive RGB split in the hk index in M22 is 
originated due to the heterogeneous heavy elemental abundances, in particular calcium. In 
our previous work, we demonstrated that the RGB stars in th e Ca-s gr oup in M22 are in 
fact Ca-rich, based on seven RGB stars studied by Brown & Wallerstein (1992). 


Here, using the e xpan ded sample of RGB stars studied by iNorris fe Freeman (119831 1 . 
Da Costa et al.l ( 20091 1. and Mar ino et ah (1201111 . we compared the differences in the elemen¬ 
tal abundances between the Ca-w and Ca-s groups in M22, providing more observational 
evidences that the two populations of stars are intrinsically different. 


4-2.1. Comparisons with Norris & Freeman (198&) 


Hesser et al.l (119771) and iHesser fe Harris (119791) first noticed that M22 have RGB stars 
with heterogeneous light elemental abundances using the DDO photometry and the low 
resolution spectroscopy, respectively, leading them t o su gg est t hat M22 may share the similar 


elemental abundance anomalies of to Gen. Later Norris fe Fr eeman (119831 1 performed a 


spectroscopic survey study for about 100 RGB stars in the cluster and they showed that the 
calcium abundance is correlated with the CN band and the G band strengths in RGB stars 
in M22. 

In Figure H31 we sh ow a plot of the calcium index A(Ca) as a function of V magnitude 
[see Equation (2) of [N orr is fe Freeman for the definition of the A(Ca) index]. In the Figure, 
we also show the least square fit to the Ca-w RGB stars, which can provide a baseline, 
and we calcul ate t he deviations of individual stars from the fitted line. As discussed by 


Norris fe Freeman (119831 ). the residuals of individual stars from this fitted line against the 


V magnitude may correct the temperature and the surface gravity effects on th e absorpt ion 
strengths to the first order. The figure shows that the A(Ca) index of INorris fe Freeman 


(119831 ) is correlated well with our hk index, which should not be a surprise because both 
quantities measure the same Ca II H and K line strength of RGB stars photometrically and 
spectroscopically. 


4-2.2. Comparisons with Da Costa et al. ( 2009 ) 


Da Costa et ah] (2009) studied M22 RGB stars using intermediate resolution spectra and 


they obtained infrared Ca II triplet strengths E(CaT) (= IL 8542 +H 8662 ) for membership RGB 
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stars, finding a substantial intrinsic metallicity spread in M22 RGB stars. The advantage 
of using the infrared Ca II triplet is that these lines arises from a lower energy level of 1.70 
eV, while the lower energy level of Ca II H and K lines is 0 eV. Therefore, the infrared 
Ca II triplet lines do not suffer from interstellar reddening contamination! Also, since the 
absorption lines from other elements around the infrared Ca II triplet lines are sparsely 
populated, contamination from absorption lines of other elements would be less severe. 

In the upper panel of Figure [GO we show the plot of E(CaT) against V — Vhb, the V 
magnitude difference from the HB. In the figure, we also show the least square fit to the 
data with a dotted line. The residuals around this fitted line will correct the temperature 
and the surface gravity effects on E(CaT) to the first order. As can be seen, the RGB stars 
in the Ca-s group have larger E(CaT) values at a given V — Rhb magnitude than those in 
the Ca-w group, consistent with our results that the Ca-s RGB stars are more calcium rich 
than the Ca-w RGB stars in M22. 


4-2.3. Comparisons with lMarino et all 120011. \2011 1) 


Finally, we make use of results from high resolution spectroscopy by Marino et al. (2 0111) . 
In Figure ITol we show plots of the [Ca/H] versus [Fe/H] and the [Na/Fe] versus [O/Fel for 
RGB stars in the Ca-w and the Ca-s groups. Similar to the results in Lee et al. u2009allbl ). 
the plot of [Ca/H] versus [Fe/H] shows that the RGB stars in the Ca-w group have lower 
calcium and iron abundances than those in the Ca-s group, confirming that our hk index 
measures the calcium abundances of RGB stars in M22. The plot of [Na/Fe] versus [O/Fe] 
shows that each group based on our hk index exhibits its own Na-0 anti-correlation. This 
was also pointed out by [Marino et ah (2011), who found two separate Na-0 anti-correlations 
in M22 RGB stars with different s-process elemental abundances. Note that the separation 
in the hk index in our current work is equivalent to that in s-proces s elemental abundan ces 
bv lMarino et al.l (120111) . These separate Na-0 anti-correlations made lMarino et al.l (120111) to 
suggest that M22 may be composed of the merger of two clusters. In fact, NGC 1851 shows 
the simil ar be havior in the Na-0 anti-correlations between the two distinct RGB populations. 
Carretta et al.l (12010) suggested that the simple explanation for the observational evidences 
of NGC 1851 would be a merger scenario, preferentially in a dwarf galaxy environment. 


Lee et al. (l2009al) showed that the differential reddening effect on the Ca II H and K lines across 


the cluster does not affect on the RGB split in M22. We demonstrated that, for example, the differential 
reddening effect of A E(B — V) « 0.3 mag results in A Ca < 0.01 mag, too small value to fully explain the 
RGB split in the hk index in M22. 
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However, it is worth emphasizing is the fact that the mean values of the Ca-w RGB 
stars appear to have a higher oxygen and a lower sodium abundances, representative of 
the primordial elemental abundance pattern, than the Ca-s RGB stars, consistent with 
the result bv iMarino et al.l 1)201 lh . Also lighter eleme nts C, N and heavy elements Ca and 
Fe abundances are different between the two groups (IMarino et al.l 20111) . Therefore, the 
feedback on the overall chemical evolution between the two groups can not be completely 
neglected. 


From comparisons of our hk photometry with three independent works by others, we 
conclude that our hk index measurements for RGB stars in M22 are very consistent with 
calcium abundance measurements by others from Ca II H and K, infrared Ca II triplet 
lines and visual Ca I lines. Therefore, the hk index can be securely and efficiently used to 
distinguish the difference in the calcium abundances at a given luminosity and furthermore 
the multiple RGB populations in GCs. 


4-2.4- On the internal 5CN vs 5CH positive correlation 


Lim et al.l (j2015[) presented low resolution spectr oscopy for M22 RG B stars and con- 


firmed the CN-CH positive correlation discovered by iNorris fc Freeman (119831 ) (see also 
Anthonv-Twarog et al. 19951) FI Here, we discuss more about the metallicity effect on the 


CN-CH positive correlation between the two RGB groups in M22. 

The variations in lighter elements in GC stars have been known for more than three 
decades, pioneere d bv ICohenl (119781). and it i s now gener ally believed to be rooted in the 
primordial origin (ICottrell fe Da Costal 11981c iKraftl Il994f) . The CN-CH anti-correlation is 
one of well-known characteristics in the Galactic GC s yste m. Ontyprobablc explanation for 
origin of this CN-CH anti-correlation can be found in Grundahl et ah (2002). They nicely 
demonstrated that there exist a CH-NH anti-correlation and a CN-NH positive correlation 


4 It is worth mentioning that iKavser et al 
M22. Examinations of their data led Pancino et al 


( 20081) did not find any CN-CH positive correlation in 


( 201 0) to attribute this to low S/N ratios of spec¬ 
tra by Kayser et al. (2008 [ and the differ ential reddening effect. However, there is other possibility that 
the measu rements by Kavser et al. ( 2008 ) may be not correct. For example the CH strength measure- 
ment s by Kayser_etalJ ( 2008 1 are unlikely very large for their ad opted definit i on of the CH strength 


from Harbeck et al 


be unreliable. 


suits by 


( 20031), indicating that the results presented by IKavser et al. (120081 ) are suspected to 


Pancino et ah ( 201 0) examined the 5CN versus the <5CH distribution using the data by 


Kavser (dahy 2008). fi nding neith er a p ositive correlation nor an anti-correlation, contradictory to the re- 


Norris fc Freeman 


(j 19831 ) or bv iLim et al.l (|2015l ). 
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among RGB stars in NGC 6752, suggesting that the nitrogen abundance is inversely corre¬ 
lated with the carbon abundance. 


The CN-CH positive correlation shown in Figure 13 of Lim ct ah (120151 ) appears to show 
two separate CN-CH anti-correlations among stars in each group at a first glance. It would 
have been very surprising if there exists a single CN-CH anti-correlation between the Ca-w 
and the Ca-s groups in M22. It is thought that a CN-CH positive correlation superposed on 
two separate CN-CH anti-correlation in M22 can be expected naturally if M22 is composed of 
two groups of stars with heterogeneous metallicities. It is because the positive correlation in 
the <5CN versus 5CH relation can occur when one com pares t he <5CN and the <5CH of different 
groups of stars with heterogeneous metallicities (see Norris et al. 1981, for the definition of 
<5CN and 5CH). 

In order to investigate the metallicity effect on the positive correlation between the CN 
and CH band strengths, we used the homogeneous measurements of the CN and CH band 
strengths for stars in 8 GCs observed in the course of Sloan Digital Sky Survey (SDSS) by 
Smolinski et ahl (2011). Since they provided magnitudes of stars in the SDSS photometric 
system , we co nv erted their SDSS photometry to the Johnson system using the relation given 
by Jester et ah ( 2005 1. 

V — g — 0.59(<7 — r) — 0.01, (2) 


and we calculated the V — Vhb using the Vhb values given by [Harris) (1199611 . 


In Figure US! we show distributions of the CN(3839) and the CH(4300) indexes against 
V— Vhb for 8 GCs. As shown, the CH(4300) strengths are well correlated with the metallicity 
of GCs at a given magnitude (for example, the RGB stars around the level of the HB mag¬ 
nitude shown with vertical dashed lines), while the CN(3839) strengths exhibit rather weak 
correlation against metallicity mainly due to the multi- modal distributions_of the CN(3839) 
strengths of individual GCs (for examp l e, see Figure 4 of lSmolinski et al.ll201 lh . Also shown 
ar e M22 and NGC 288 fro m iLim et ahl (1201511 , whose measurements do not agree with those 
of ISmolinski et ahl (120 111) due to heterogeneous definition of the band strengths between 
the two wor ks. Note tha t ISmolinski et ahl (1201111 used the definiti on of the CN(3839) b and 
strength by Norris et ah ( 1981 1 for RGB and SGB stars, that by Harbeck et ah ( 2003h for 


MS st ars and the definition of the CH(4300) band s trength by iLeel (119991) , while ILim et ah 
(2015) used the definitions by Harbeck et ah (120031 ) for both bands. As a consequence of 


heterogeneous definition of the m olecule band stren gths, the extent of the CN(3839) band 


stren gths for M22 and NGC 288 bv lLim et ahl (120151) is great er than tho s e bv ISmolinski et ah 
( 2011 1. Also, the CH(4300) values for M22 and NGC 288 by him et ah ( 20151) are > 1.0 and 


they are out of range in the figure. Therefore, a direct comparison between the two works is 
not feasible. 
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In the figure, the thick solid lines are the common lower envelopes for 8 GCs. Due to 
the difference in the mean metallicity among GCs, these common lower envelopes may not 
correct both the temperature and the surface gravity effects simultaneously on the CN(3839) 
and the CH(4300) band strengths for 8 GCs. This is exactly what expected to happen when 
one compares two groups of RGB stars in M22 with different mean metallicities using a 
single fitted line, since the two groups of RGB stars in M22 have very different metallici- 
ties. In Figure [T71 we show the distributions of the hCN and the 5CH against metallicity, 
where the 5CN and the £CH are defined to be the differences in the CN and CH between 
individual stars and the common lower envelopes. As expected from Figure [161 the 5CN 
and the 5CH have substantial gradients against metallicity, 0.169 ± 0.031 and 0.077 ± 0.010 
mag/dex, respectively. It is believed that these positive slopes are the source of the CN-CH 
positive correlation when one compares the <5CN and the <5CH of stars with heterogeneous 
metallicities. 


In Figure HU we show the CN(3839) and the CH(4300) band strengths by Lim ct ah 


( 20151 ) against V magnitude. Also shown are the lower envelops for the Ca-w (the blue 
dotted line) and the Ca-s (the red dotted line) groups. Stars with —1.0 < V — Vhb < 0.5 
mag are shown with plus signs (Ca-w) and crosses (Ca-s). Note that the lower envelope for 
the Ca-w group can be served as the common lower envelope for both groups. 


Figure [191(a) shows a distribution of the <5CN versus the <5CH of M22 RGB stars, where 
the <5CN and the <5CH are calculated with respect to the common lower envelope (i.e. the blue 


dotted lines in Figure [IH]) . This plot is similar to Figu re 13 of 


CN-CH positive correlation can be seen. Noted that Lim ct ah 


Lim c 


al.l (120151) . where the 


(120151) used the mean fitted 


line, not the lower envelope used in this analys is or those of others (for example iNorris et al. 


198ll ; iKavser et al.ll2008l: ISmolinski et al.ll201 ll ) and detailed distributions of individuals stars 


are slightly different from those presented here. 


In Figure [191(b) and (c), we show distributions of the <5CN versus <5CH for the Ca-w 
and the Ca-s groups using their own lower envelopes. Each plot appears to show its own 
CN-CH anti-correlation, a well-known characteristics of normal GCs. This is also consistent 
with the prese nce o f the separate C-N anti-correlations in each of the two s-process groups 
(Marino et ah 2011). Since its own lower envelope for each group can correct the metallicity 
effects, the 5CN and the <5CH values calculated from its own lower envelope reflect the relative 
CN and CH band strengths at a given metallicity and surface gravity. As a consequence, the 
CN-CH positive correlation with an apparently ambiguous origin disappears. In this regards, 
the most simplest explanation would be a merger of two GCs with different metallicities. 
Each group of stars could have been chemically enriched independently and exhibits its own 
CN-CH anti-correlation as observed in normal GCs in our Galaxy. 
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Finally, it is worth point i ng ou t again that the extents in the hCN and the <5CH of M22 
using the data by Lim ct ah (120151 ) are greater than those can be expected from the slopes 
in the metallicity gradients as shown in Figure [T7J This is partially due to the different 
definitions in the CN and the CH band strengths. It would be highly desirable to re-analyze 
spectra in a homogeneous manner in the future. 


4-2.5. Spread in nitrogen abundances 


Yong et ah (12008 ) investigated nitrogen abundances of RGB stars in NGC 6752 using 
NH molecule lines. They noted that Stromgren u passband contains the NH molecule band 
at A 3360 A and they showed that the c\ index are correlated with the nitrogen abundances 
of individual stars. They devised an arbitrary empirical index cy [= Ci — (b — y)} and derived 
a well correlated relation between the nitrogen abundance and the cy index of RGB stars at 
about Vh b °f NGC 6752, [N/Fe] oc 16.11 x cy. 


In Figure 1201 we show cy versus V CMDs of M22 RGB stars, showing broad spreads in 
the cy index in both groups. Assuming Gaussian distributions, we calculated the standard 
deviation around the mean cy values and obtained a(cy) ~ 0.03 mag at Vh b for both groups 
of stars in M22. Since the metallicity of M22 is comparable to that of NGC 6752 and if the 
spread in the NH absorption strengths in individual RGB stars is solely responsible for the 
spreads in the cy index, cr(cy) ~ 0.03 mag corresponds to a [N/Fe] fa 0.5 dex in each group 
of stars. 


Note that our photometric estimate of the spread in the nit rog en _abundances in M22 
RGB stars appear to be slightly larger than those of Marino et ah (201l|), who obtained 
cr[N/Fe] ~ 0.2 dex for their s-poor and s-rich groups based on the observations of 14 stars. 


One of the remarkable features can be seen in Figure [20] is that the RGB bump magni¬ 
tudes of both populations appear to be different. In the right panel of Figure [201 we show 
the cumulative luminosity functions (LFs) for the Ca-w and the Ca-s populations. We show 
the V magnitude levels of the RGB bump, V , ump , where the slope of the luminosity function 
changes abruptly, for each population by arrows. We obtained V\, ump = 13.91 mag for the 
Ca-w population while 14.06 mag for the Ca-s population. 


4-2.6. RGB bump magnitude: Vh ump versus A hk 

Here, we explore more details on the RGB bump in M22. As already shown in the left 
panel of Figure [7] it can be seen that the V magnitude level of the RGB bump appears 
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to be positively correlated with the A hk values of RGB stars at around V ~ 14.0 mag. 
Especially, due to its high stellar number density, the tilted 14 urnp magnitude level against 
A hk for the Ca-w group is conspicuous, while the tilted HB structure in the Ca-w group 
cannot be seen in the cy versus V CMD as shown in Figure [201 It should be worth noting 
that neither this tilted V magnitude against A hk in the Ca-w group nor the difference in 
the bump magnitu de betw een the two groups can be explained by the differential reddening 
effect (see also, Lee et aL 2009al ). 


The difference in the mean A hk values between the Ca-w stars with A hk < 0 mag and 
those with A hk > 0 mag at V — 13.95 mag is A hk = 0.04 mag. Again, assuming that 
the conventional interstellar reddening law by [Anthony-Twarog et alj (19911) is applicable, 
E(hk) = —0.16 xE(b — y) ~ —0.12 xE(B — V), and that the differential reddening effect 
across the cluster is the main reason for showing a broad RGB sequence in the hk index, this 
0.04 mag difference in the A hk index can be translated into Ay = —1.05 mag, in the sense 
that the Ca-w RGB stars with A hk > 0 mag should be about 1 mag brighter than those 
with A hk < 0 mag. This is the opposite of what observed in M22. As shown in Figure I2T1 
the V magnitude for the RGB bump of the Ca-w stars with A hk > 0 mag is about 0.06 
mag fainter than that of the Ca-w stars with A hk < 0 mag (see also Figure 1201). Also, if 
this 0.04 mag difference in A hk is originated solely by the differential reddening across the 
cluster, not by the difference in the chemical compositions, the (b — y ) color of the Ca-w 
group should have been broadened by m 0.25 mag. However the RGB sequence in M22 is 
not as broad as 0.25 mag as shown in Figure Of Therefore, the difference in the RGB bump 
magnitude is not mainly due to differential reddening effect but due to the difference in 
chemical compositions. 

It is well known that at a given age the RGB bump becomes fainter with increasing 
metallicity and with decreasing helium abunda nce, due to changes in the envelope radia¬ 
tive opacity (see for example, C assis i & Safaris 20131) . The fainter Vbump magnitude with 
increasing A hk can be interpreted that A hk increases with either increasing metallicity or 
decreasing helium abundance. One can quantitatively estimate the effect of metallicity and 
the h elium abundance on the RGB bump luminosity flBiork fe Chaboverll2006l: IValcarce et al. 


2012 ). 


In their Table 2, Biork fe Chabover ( 20061 ) presented the absolute V magnitude of the 
RGB bump with metallicity and we obtained AMyb U mp/A[Fe/H] « 0.93 m ag/dex for 13 Gyr. 
The effect of helium abundances on the RGB bump can be found in Valcarce et ah| (2012). 
From their Figure 9, we obtained Amboi ~ 2.5 x AT' for the isochrones with Z — 1.6 x 10~ 3 
and 12.5 Gyr. If the 0.06 mag spread in the RGB bump magnitude in the Ca-w group is 
solely due to metallicity, one expects to see the spread in metallicity by A[Fe/H] ~ 0.06 
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dex. Similarly the 0.15 mag difference between the Ca-w and the C a-s gro ups can translate 
into A[Fe/H] « 0.16 dex that agrees well with what Marino et al.j (2011j) obtained for the 
s-poor and the s-ricli groups, A[Fe/H] ~ 0.15 ± 0.02 dex. However, if we take int o ac count 
the enhanced helium abundance for the Ca-s group by AY & 0.09 ( Joo & Leel 120131) . the 
RGB bump magnitude for the Ca-s group should be ~ 0.08 (= 0.23 — 0.15) mag brighter 
than that for the Ca-w group, which is not the case for our observation. The formation 
of the EBHB stars, as progenies of the Ca-s RGB stars, may require the enhanced helium 
abundance but the difference in the RGB bump magnitude is difficult to explain in this 
scenario. Perhaps, may this suggest that only a fraction of the Ca-s RGB stars have been 
evolved into the EBHB? 


4.3. Centers 


With the bright RGB stars selected above, we measured the centers of each popula¬ 
tions using three different methods; the arithmetic mean, the half-sphere and the pie-wedge 
methods. 


Using the coordinate of the center of the cluster measured by Goldsbury et al. (120101 ) 
as an initial value (a = 18:36:23.94 and 6 = —23:54:17.1; J2000) we chose RGB stars in 
each population within 300 arcsec from the center of the cluster and we calcul ated the mean 
values for each group. The half-light radius of M22 is about 200 arcsec (IHarria 19961) and the 
circle with a radius of 300 arcsec can contain more than 70% of RGB stars in our sample. R 
is believed that using this radius is sufficient enough to derive the coordinate of the center 
of each group with decent accuracy. 

We obtained the offset values with respect to the coordinate of center of the cluster by 
Goldsbury et ah . (ARA, A Dec) = (3'.'8, 0'.'4) for the Ca-w group and (0"8, 5'.'6) for the Ca-s 
group, resulted in an angular separation of 6 arcsec between t he ce nters of both populations. 
Note that the core radius of the cluster is about 80 arcsec ( Harris! 19961) . and the angular 
separation between the centers of the two populations is relatively small. 


Similar to the simple mean calculation, we use the coordinate of the cluster by Goldsbur y et al, 
as an initial value, we chose RGB stars in each population within 300 arcsec from the center 
of the cluster. Then we divided the sphere into two halves by assuming the radial symmetry 
in the distribution of RGB stars in M22. We compared the number of RGB stars between the 
two halves by rotating the position angle by 10 degree at a fixed coordinate of the center and 
we obtained the differences in the number of RGB stars between both halves. We repeated 
this calculation with varying coordinates of the center and we derived the coordinates of the 
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centers of each group with the minimum difference in the number of RGB stars between the 
two halves. We found (775, 17 6) for the Ca-w group and (—174, —174) for the Ca-s group, 
slightly different from those from the simple mean method. The angular separation between 
the two groups is about 9 arcsec and, again, it is relatively small compared the core radius 
of the cluster. 


Finally, we applied the pie-wedge method. Using the coordinate by Goldsbu rv et ah 


( ,20101) . we divided the sphere of a radius of 300 arcsec into 12 different slices. Then we 
compared the number of stars in the opposing distribution. We repeated this calculation 
with varying coordinates of the center and obtained the center of each population with the 
minimum differences. We found (7'!7, —2 , . / 8) for the Ca-w group and (—4"9, 675) for the 
Ca-s group. Again, the angular separation between the two groups is about 16 arcsec and it 
is still relatively small compared with the core (~ 80 arcsec) or the half-light (~ 200 arcsec) 
radii of the cluster. 


We show our results in Tableland Figure [22] We conclude that the coordinates of the 
centers of each population are slightly different. However the differences in the coordinates 
of the center from various methods do not appear to be substantially large to claim that the 
centers of the RGB distributions of two populations are distinctively different. 


4.4. Spatial Distributions 

4-4-1- Red, giant branch stars 


In Figure [23l we show the number ratios of RGB stars in the Ca-s group to those in the 
Ca- w gr oup against the radial distance from the center using the coordinate of the center 
by Goldsbu rv et ah 020101 1. It is worth noting that there appears to exist a weak gradient of 
the number ratio between the two populations against the radial distance from the center, 
indicating that the distribution of the Ca-w group is slightly more centrally concentrated. 
However, the difference in the mean number ratio are in agreement within the measurement 
errors from the central part of the cluster up to ~ 800 arcsec from the center, about four 


times of the half-light radius of M22, similar to that can be seen in NGC 1851. Milone et ah 


f 20080 showed that the number ratio between the two distinctive populations in NGC 1851 
does not appear to vary up to 1.7 arcmin from the center, which is about 2 times of the 
half-light radius of the cluster]^ The flat number ratio between the multiple populations 


Carrett a et all (120111 ) claimed that the metal-poor RGB population in NGC 1851 is more centrally 


concentrated than the the metal-poor population from the incomplete sample of their spectroscopic study 
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in M22 and NGC 1851 is in sharp contrast with that in uj Cen. Bellini ct ah (2009) found 
that both intermediate-metallicity and the metal-rich RGB populations are more central ly 
concentrated than the metal-poor RGB population in u Cen (see also Sollima et ah 20071 ). 


The relaxation time at the half-mass radius for uj Cen is about 10 Gyr ([Harris 1990 ) and 
it is substantially large that the later generation of stars in u Cen may not have sufficient 
time to be homogenized completely. As a consequence, the later generation of stars fails 
to achieve similar radial distribution as the early generation of stars in the cluster. On the 
other hand, the relaxation times at the half-mass radius for M22 and NGC 1851 are about 
1.7 Gyr and 0.7 Gyr, respectively, significantly smaller than that of u Cen. If the second 
generation of stars in M22 and NGC 1851 formed from gaseous ejecta expelled from the first 
generation of stars, they might have enough time to become relaxed systems. In this regard, 


Decressin et ah (12008!) claimed that stars initially centrally concentrated toward the center 


need about the two relaxation times to achieve a completely homogenization throughout the 
cluster and any radial difference among multiple stellar populations would be erased over 
the dynamical history of old GCs@ However, it is also possible that the merger of two GCs 
can maintain a flat RGB number ratio in M22. 


We show a comparison of cumulative distributions of RGB stars in both groups in the 
lower panel of Figure |23j As can be seen, the radial distribution of RGB stars in the Ca-w 
group is slightly more centrally concentrated. We performed a Kolmogorov-Smirnov (K-S) 
test and we found the probability of being drawn from identical populations is 0.1%, with a 
K-S discrepancy of 0.10, indicating that they have different parent populations. As we have 
demonstrated earlier, the difference in the radial distribution is not likely due to the held 
star contamination. 


Although small, the difference in the distributions of RGB stars is difficult to understand 
within the current theoretical framework of the GC formation. If the second generation of 
stars formed out_of gas expelled from the first generation of stars in the central part of the 
cluster (iBekkii 120101 . see for example), the second generation of stars should have a more 
centrally concentrated distribution in the early history of the cluster. Also importantly, if 
the complete homogenization had been achieved within a couple of relaxation times, the two 
different populations in M22 are expected to have indistinguishable radial distributions. 


with larger radial distance from the center, « 12 aremin. 


6 The work of Vesperin i et al. ( 20131 ) is also worth mentioning. They performed IV-body numerical 
simulations and showed that the time required to have a flat number ratio between the first and the second 
generations can be larger than 45 half-mass relaxation time for their rlO system, where rlO refers the ratio 
between the two half-mass radii (see their Figure 7). 
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We speculate that the difference in the radial distributions of RGB stars may suggest 
that either (or perhaps both) of followings is correct; (i) The Ca-s population did not form 
out of the gaseous ejecta of the first generation of stars. Or, (ii) the complete homogenization 
may require much longer time scale than that proposed by Decressin et al. (2008). In this 
regard, the merger of two GCs can naturally explain the observed radial distributions of 
RGB stars. 


We turn our attention to the spatial distributions of RGB stars in both groups. In the 
top panel of Figure [2H we show the projected spatial distributions of RGB stars in the Ca-w 
and the Ca-s groups in M22. In the Figure, the offset values of the projected right ascension 
and the declination i n the units of arcsec were calculated using the transformation relations 
in 


van de Yen et ah 020061 ). 


. „ . 648000 . . . 

AKA = -coso sm Aa, 

7T 

648000 

A Dec = -(cos o cos oo — coso sm o 0 cos Act), 

7r 


(3) 


where Aa = a —a 0 and A<5 = 5 — 5 0 , and a 0 and 5 0 are the coordinate of th e center. For our 


calculations, we adopted the coordinate for the cluster center measured by Goldsburv et ah 

fhoioh . 


In the lower panel of the Figure, we show smoothed density distributions of RGB stars 
along with iso-density contours for each population. For the smoothed density distribution 
of each population , we appl ied a fixed Gaussian kernel estimator algorithm with a FWHM 
of 70 arcsec (1 Silverman! 1 9861 ). We show the FWHM of our Gaussian kernel in the lower left 
panel of the figure. To d erive the iso-density contour for each p opulation, we applied the 
second moment analysis (iDodd fc Ma,cGillivarvlll986l; IStonel Il989l) . In the Figure, we show 
the iso-density contour lines for 90, 70, 50, and 30% of the peak values for both populations. 


At a glance of the figure, the distribution of the Ca-s group appears to be more spatially 
elongated than that of the Ca-w group does. In Table |6l we show the axial ratio, 5/a, and 
the cllipticity, e (= 1 — 5/a), and we show radial distributions of the axial ratio and the 
cllipticity of the Ca-w and the Ca-s groups in Figure 1251 As shown, the radial distributions 
for both groups agree rather well up to ~ 130 arcsec (~ 1.5 core radii of the cluster) from 
the center but two distributions bifurcate at the larger radial distance, in the sense that the 
Ca-s group is more elongated. 
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4-4-2. Horizontal branch stars 


In recent years, the nature of the multiple stellar populations of the HB stars in GCs 
has been emerged th r ough the deta i led sp ectroscopic studies (e.g. iMarino et al.ll2013l. l2014al: 


Gratton et al.ll2014J) . IMarino et al.l ( 2013 1 studied 7 HB stars in M22 and they found that all 


their HB stars are Ba-poor and Na-poor, equivalent to the s-poor RGB population, leading 
them to suggest that th e pos ition of a HB star is strictly related to the chemical composition. 
More recently, Gratton et ah (2014) investigated the chemical compositions of more than 90 
HB stars in M22, using the multi-object spectroscopy facility. Their results showed that 
there appears to exist 3 different groups of HB stars in M22. However, they did not find 
any severely O-depleted HB stars. They claimed that the progeny of the severely O-depleted 
RGB stars may evolve into the hotter part of the HB, where the surface temperature of these 
HB stars are beyond the temperature range of their sample. 


In addition to metallicity and age, it is generally b elieved that the HB morphology 


of GCs can be goyerned by the helium abundances (e.g. iD’Antona et al.l l2002[ |Joo fe Lee 


20131; iMilone et al.l 1201411 and there als o exist direct spectroscopic observational evidences of 
helium enhancement in GC HB stars (IMarino et al.l l2014al: iGratton et al][2014j) . 


In their recent work, Joo Lee (120131 ) investigated M22 HB populations using the 
synthetic population models. Although simplistic, they divided M22 HB populations into 
two groups, the BHB and the EBHB, and they delineated a connection between the RGB 
and the HB of M22 (see their Figure 13). Here, we explore more details on the number ratio 
and the spatial distribution of the HB populations in M22 using our wide field data. Note 
that the hk index is not an useful tool to study the BHB or the EBHB stars, since the Ca II 
H and K lines become suppressed and the absorption strength of H e at A 3970 A reaches its 
maximum at A0 and can contaminate the adjacent Ca II H line at the high effective surface 
temperature of BHB or EBHB stars. 


In Figure [23 we shows CMDs of the HB region of M22. As for the RGB populations, 
we used all available multi-color information in order to remove the contamination from the 
off-cluster field star population in our study. Also shown in the figure is the cumulative LF 
of the HB stars, where the slope in the cumulative LF of HB stars becomes briefly zero at the 


location of the 1 

IB ga 

p at V = 15.97 mag. This magnitude is about the same magnitude level 

that 

Joo & Lee 

(2013 

) separated the BHB and the EBHB populations in M22. Furthermore, 

Joo & Led ( 

20R 

1) claimed that the BHB stars are the first generation of stars (corresponding 


to the progeny of the Ca-w RGB population in our notation), while the EBHB stars the 
second generation of stars (the Ca-s population) of the cluster. 


In our photometric data, the observed number ratio between the BHB and the EBHB 
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stars becomes n(BHB):n(EBHB) = 80:20 (± 9), marginally in agreement with those of the 
RGB populations, n(Ca-w):n(Ca-s ) = 69:31 (± 3), and the SGB populations, n(bSGB):n(fSGB) 
= 69:31 (± 6), without the correction for the evolutionary effects. With the correction of the 
evolutionary effect, the HB number ratio does not agree with that of SGB, n(bSGB):n(fSGB) 

= 62:38 (± 6). Note that our ffB number ratio includes the completeness correction that 
shown in Figured We calculated the average completeness fraction as a function of V mag¬ 
nitude using those from the central and outer part of the cluster as shown in Figure [3J. and 
we applied the correction factor to individual HB stars by interpolating the average com¬ 
pleteness fraction against V magnitude. As discussed above, our photometry is complete 
down to ~ 18.5 mag. Therefore, only the the EBHB population is insignificantly affected 
(~ 0.6 %) by the completeness correction. 

In the upper panel of Figure [23 we show the number ratios of the EBHB to the BHB 
populations against the radial distance from the center. Although, the HB number ratios 
against the radial distance are in agreement within measurement errors up to ~ 3 times 
of the half-light radius of the cluster, the mean HB number ratio appear to have a weak 
gradient against the radial distance in the sense that the BHB stars are slightly more centrally 
concentrated than the EBHB stars, similar that can be seen in the RGB stars in Figure [221 
In the bottom panel of the figure, we show the cumulative distributions of the BHB and 
the EBHB populations against the distance from the center. The BHB population is more 
centrally concentrated than the EBHB population, qualitatively consistent with that can be 
seen in the RGB populations. The K-S test indicates a probability of 1.8% that two HB 
populations are drawn from the same parent population, suggesting that they have different 
parent populations as can be seen in the RGB stars. 


As have done for the RGB and SGB populations, due to different masses and chemical 
compositions, different lifetimes of individual HB stars should be taken into account. We 
comp ared t he lifetimes between the BHB and the EBHB populations using the HB tracks 
by Joo fe Leel ( 2013 ) and obtained the correction factor in terms of the number of HB stars, 
n(BHB):n(EBHB) = 56.5:43.5 (± 4.0). Because the EBHB stars live longer, one naturally 
expects to observe a larger number of EBHB stars if initially even numbers of stars are 
distributed both on the BHB and on the EBHB locations. 


In addition, the number of HB stars evolved from the Ca-w and the Ca-s RGB pop¬ 
ulations will be slightly different since the RGB-tip masses between the two groups will be 
different. If the Ca-s population has the higher helium abundance by AY = 0.09 (Joo fe Lee 


2013 1). the evolution of stars in the Ca-s population would be faster than those in the Ca-w 


population due to increasing mean molecular weight, resulted in s maller RGB-tip masses 
in the Ca-s population. With the HB mass dispersion adopted by Joo fe Lee 12013 ). om 
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= 0.023M q for both populations, we calculated the number of stars within 4 xctm from the 
RGB-tip massed of both isochrones using the same evolutionary population synthesis model 
that we constructed previously. For this purpose, we populated 10' artificial stars in each 
group using the Salpeter’s IMF and we generated 50 different sets for both groups. We 
obtained the correction factor, n(BHB):n(EBHB) = 58.9:41.1 (± 0.1), in the sense that the 
observed number of HB stars that evolved from the Ca-s population gets larger since the 
RGB-tip mass of the Ca-s is smaller owing to its faster evolution from the MS through the 
RGB phases. If we take these effects into consideration, the number ratio between the BHB 
and the EBHB stars becomes n(BHB):n(EBHB) = 88.2:11.8 (± 9.8), significantly different 
from those of the RGB and the SGB populations. It should be noted that this result is based 
on the assumption that the BHB stars are the progeny of the Ca-w RGB sta rs w hile the 
EBHB stars are the progeny of the Ca-s RGB stars as proposed by (loo & Lee (120131 ). 


Our result may sug gest t hat the connection between the RGB and HB populations may 
not be as clear as that Joo fe Lee (120131 1 claimed. Note that IJoo fe Lee (120131 ) adopted the 
number ratio between the BHB and EBHB of 70:30 in their synthetic HB model construction 
and they claimed that their HB number ratio does not disagree with the previous estimates 
of RGB or SGB number ratios between the two populations by others. It is thought that 


Joo & Lee (2013 ) overest imated the EBHB fraction in their calculation compared to our 


observation. As Joo & Lee (j2013i ) suggested, if the EBHB population is the only progeny 
of the second generation of the RGB population (i.e. the Ca-s population), the discrepancy 
in the number ratios between the two populations in the HB and SGB/RGB is difficult to 
explain. For instance, it is evident that the R value, the number ratio of the HB to the RGB 
stars (jCa puto et al.l 1 9871 ). of each population is against the theoretical expectation. The 
observed R value of a simple stellar population is the measure of the helium abundance in 
the sense that the R value increases with the helium abundance. We define the mo diked R 
value which includes RGB stars from Vp b ( = 14.15 mag for M22) to 12.0 mag (not the tip 
of the RGB sequence), since we already selected bright RGB stars with 12 mag < V < 16.5 
mag in both populations as discussed above. Then we have the modified R values of 1.73 
± 0.14 for for the BHB versus the Ca-w RGB populations and 1.21 ± 0.17 for the EBHB 
versus the Ca-s RGB populations. If the EBHB and the Ca-s RGB populations are the 
second generation of the stars in the cluster and they formed out of interstellar gas enriched 
in helium by the first generation of stars, one would naturally expect to have a larger R 
value for the EBHB and the Ca-s RGB populations. However, the observed R value for the 
EBHB and the Ca-s RGB populations is smaller than that for BHB and the Ca-w RGB 


7 This mass range is rather arbitrary, but the final result does not significantly depend on the adopted 
mass range if the same mass dispersions for both groups are maintained. 
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populations and the observe R values are opposite of the theoretical expectation. Our result 
may suggest that, like other lighter elemental abundances, there may exist a spread in the 
helium abundance of the second generation of stars and at least about the half of the Ca-s 
RGB stars must have evolved into the BHB sequence (which requires no or little helium 
enrichment) and a significant fraction of the Ca-s RGB s tars failed t o bec ome EBHB stars 
(which requires the helium enrichment of AY « 0.09, see Joo & Le^ 120 130 in M22. 


ft would be very interesting to point out the recent result bv lGratton et ahj (2014). They 
carried out the spectroscopic survey of the BHB stars in M22 and found a significant fraction 


(~ 18%) of the m etal-rich helium-enhanced BHB stars in their sa mple (Group 3 c 


bv lGratton et all) . If we consider that the Group 3 BHB stars bv iGratton et al.l (120141) are 


esignated 


the progeny of the Ca-s RGB stars, the number ratio between the second generation and 
the first generation of HB stars becomes 34:66 and this number ratio is in good agreement 
with those from the RGB and the SGB stars. Perhaps, this is also hinted by the separate 
Na-0 distribution of M22 RGB stars, indicative of internal He spread in both the Ca-w and 
the Ca-s groups. If so, the Ca-s population of M22 may be exactly what we observe in 
NGC 6752, which contains both the BHB and the EBHB populations maintaining a rather 
monotonous metallicity distribution. 


4.5. Radial velocity and velocity dispersion 


Recently, Lane et ahj (1200911 performed a wide-held spectroscopic survey of bright RGB 
stars in M22 using AAOmega mounted on the 3.9-m Anglo-Australian Telescope. Their 
primary concern was to test Newtonian gravity by investigating velocity profiles of several 
globular clusters, including M22, at larger radii. We made use of their radial velocity data, 
which were kindly provided by Dr. Kiss, in order to investigate the differences in kinematics 
of the multiple stellar populations in M22. 


We compared our coordinates of Ca-w and Ca-s RGB stars with those of Lane et ah 
(2009) and we selected common stars matched within a radius of 2". During onr source 


matching process, we excluded s tars having uncertainties in the radial velocity measurement 
larger than 10 km s -1 tagged by [Lane, et ah ( 20091 ). Through this process, we have 208 Ca-w 
and 87 Ca-s RGB stars. 


We show the radial velocity of individual stars in each group against the radial distance 
from the center in the upper panel of Figure [28j We obtained the mean radial velocity of 
the cluster, (v r ) = —144.89 ± 6.76 km s -1 . The distributions of the radial velocity between 
the two populations are in good agreement. Also the distribution of the radial velocity of 
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each population does not vary with the radial distance from the center. 


Using these stars, we calculate th e average radial velocities and velocity dispersions for 
Ca-w and Ca-s groups as given by Pryor & Mevlanl (119931 ). assuming that each velocity n, 
is drawn from the normal distribution, 


f(Vi) = 


= exp - 7 —t -—yy , 

yj 2 ^c + (jU) l (^c+°e,i)J 


(4) 


where a £) i is the known measurement uncertainty of ty and v r , a c are the average radial 
velocity and the intrinsic velocity dispersion of each group. For our calculations, we used 
a FORTRAN program kindly provided by Prof. Pryor, obtaining v r = —144.59 ± 0.40 km 


,-i 


s " and <7 V = 4.71 ± 0.35 km s _1 for the Ca-w group and v r = —145.77 ± 0.89 km s -1 and 
<j v = 7.55 ± 0.66 km s _1 for the Ca-s group. Our results are shown in the bottom panel 
of Figure [28l The mean velocities of both groups are in good agreement, while the velocity 
dispersion for the Ca-s group appears to be slightly larger than that of the Ca-w group. 
Note that Marino et al. (l2014bi ) studied the velocity dispersions between the two SGB in 
NGC 1851 and they did not find any difference. 


4.6. Projected Rotations 


Peterson & Cudworth (1199411 first found a clear evidence of the rotation of M22 from 


both proper motions and radial velocities of bright stars in the cluster, with the amplitude 
of the mean rotation of about from ss 6 km s -1 for 1' < r < 3' to ~ 3 km s^ 1 for 3' < r < 7'. 
Their result is an unique example of obtaining the rotational velocity of GCs using proper 
motion data. Later lLane et al. (120091 . 120101 ) confirmed that M22 has a substantial rotation 
among other GCs that they studied. In this section, we investigate the mean rotation of 
RGB stars of the double populations in M22 using three methods. 


4-6.1. Radial velocity versus position angle 


First, we estimated the mean rotation of the cluster by using the radial velocity dis¬ 


tribution against the position angle of individual stars flPeterson fe Cudworth 19941 ). In 


Figure [29l we show radial velocities of individual stars in each group against their position 
angles measured from North. We chose the RGB star^l with the proper motion member- 


6 Peterson fc Cudworth (1199411 also studied the HB stars in the cluster, of which we do not make use in 
this study. Therefore, the mean velocity of the cluster and the number of stars being used in this study are 
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ship probability larger than 90% from Peterson & Cudworth (1994) and we merge their 
radial_yelocity measurements with our photometric data. We show radial velocities from 


Peterson & Cudworth (1994J) against the position angle for individual RGB stars in Fig¬ 


ure 1291(a). We performed a sinusoidal fit to the data and we obtained the amplitude of 
the mean rotation of all RGB stars is about 4.2 km s' 1 which is comparable to that from 


Peterson & Cudworth 

(1994 

), 3.8 km s 1 ( 

Mevlan & Heggie 

1997) 


Next, we investigate the mean rotational veloci t y of the Ca-w and the Ca-s groups. 
We divided RGB stars from P eterson & C udworth (Il994f ) into two groups based on the 
hk index of individual stars in Figure |T| and we show plots of radial velocities against the 
position angle of RGB stars in the Ca-w and the Ca-s groups in Figure [291(b) and -(c). 
The amplitude of the sinusoidal fit to the Ca-w group is about 3.6 km s _1 , while that to the 
Ca-s group is about 5.6 km s _1 . The amplitude of the projected rotation of the Ca-s group 
is larger than that of the Ca-w group. 



amplitude of mean rotation for the Ca-w group is about 1.5 km s _1 , while that for the Ca-s 
group is about 3.6 km s -1 . Again, the amplitude of the mean rotation for the Ca-s group is 
larger than that for the Ca-s group. 


It should be noted that the diffe r ence i n the mean radial vel ocities of the clus ter using 
the data from IPeterson fc Cudworth ( 1994 ). —149.1 km s” 1 , and lLane et al.l ( 20091 1 . —144.6 
km s^ 1 , may suggest that two sets of radial velocity data are heterogeneous. Therefore, we 
did not attempt to merge two sets of data together to investi gate t he rotati on of the cluster. 
The a mplitude of the mean rotation using the data by lLane et al.l (120091) is smaller than 
that by IPeterson fe Cudworth (119941) . Since the number o f RGB stars being used is smaller 
for the case of using data by Pete rson fe C ud worth (119941 ). it is suspected that a few stars 
with large deviations from the mean velocity can affect the results. In the upper panel of 
Figure [291 the radial velocities of three RGB stars deviate more than 3 a level from the mean 
value; IV-20 (—166.78 km s -1 ), 2-73 (—127.65 km s _1 ), and 1-27 (—127.51 km s -1 ). If we 
exclude these three stars in our calculations, the amplitudes of the mean rotations become 
2.7 km s -1 , 2.5 km s _1 and 3.3 km s -1 for all RGB stars, the Ca-w and the Ca-s group s, 
respectively. These values become comparable to those using data from L an e et al. (12009 ). 


different from those in Peterson & Cudworth ( 19941) . 
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4-6.2. Radial velocity difference between two hemispheres 


Next, we estim ated the amplitude of the mean rotation of the cluster using the method 
described by Lane et ah (120091) . Assuming an isothermal rotation, the mean rotation can be 
measured by dividing the cluster in half at a given position angle and calculating differences 
between the average velocities in the two halves. We repeated this calculation by increasing 
the position angle of the boundary of the two halves by 10°. Then the net rotation velocity 
is the half of the amplitude of the sinusoidal function in the differences between the average 
velocities in the two halves. 


We show the differences in the mean radial velocities as a function of the position angle 
(East = 0° and North = 90°) along with the best-fitting sine function in Figure [30] Our 
result for all RGB stars in the cluster is shown in the top panel of the figure. We obtained 
the mean rotation velocity of 1.4 ± 0.3 km s -1 and the position angle of the equator (i.e. 
the axis perpendicular to the rotation axis) of the r otati on of 81° (i. e. ~ North), which are 
very consistent with those of previous studies by Lane et al. J2009 . 2010), 1.5 ± 0.4 km s 
and approximately North-South, respectively. 


-l 


In the middle and the bottom panels of Figure [30] we show plots of differences in the 
mean radial velocities for the Ca-w and the Ca-s groups. We obtained the mean rotation 
of 1.0 ± 0.4 km s -1 with the position angle of 86° for the Ca-w group and 2.5 ± 0.8 km 
s^ 1 with the position angle of 73° for the Ca-s group. Again, the amplitude of the mean 
rotation of the Ca-s group is slightly larger than that of the Ca-w group is. 


Note that the mean rotations of this method are smaller than those of the sinusoidal 
fit to individual stars on the radial velocity versus position angle. As we have demonstrated 
above, the sinusoidal fit to the individual stars suspects to be vulnerable to stars with large 
velocity deviations from the mean value. This may explain the difference in the amplitudes 
of the mean rotation between the two methods. 


4-6.3. Radial velocity profile 

Finally, we examine the rotation of the cluster by using the radial velocity profile along 
the axis perpendicular to the mean rotation. In the top panel of Figure ETJ we show the 
distributions of RGB stars in the Ca-w and Ca-s groups with the radial velocity measure¬ 
ments. In the figure, blue and red colors denote blue-shift and red-shift from the mean 
radial velocity of the cluster, respectively, and the size of each circle represents the velocity 
deviation from the mean value of the cluster. We also show grids for calculating the group 
velocity and the velocity dispersion with rectangles with the width of 800 arcsec and the 
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height of 100 arcsec. We adopt the tilted angle of the grid rectangles that we obtained in 
1 14.6.21 89° for the Ca-w group and 64° for the Ca-s group. 

We calculated the the mean radial velocity and the velocity dispersion in each grid 
rectangle. Our results are shown in the bottom panels of Figure [31] Inside the half-radius 
of the cluster, the velocity dispersion of the Ca-s population appears to be slightly larger, 
consistent with our previous result presented in 1 14.51 Also, the radial velocity profile against 
the projected distance perpendicular to the rotation axis are different. 


We performed a least square fit to the line of sight velocity versus the projected distance 
perpendicular to the rotation axis. We found that the maximum velocity of 1.0 km s -1 for 
the Ca-w group and 4.3 km s -1 for the Ca-s group at the half-light radius of the cluster. 
Again, the RGB stars in the Ca-s group appear to rotate faster than those in the Ca-w 
group do, consistent with our previous results. 


We examined the anisotropy parameter of the cluster (see, for example, Bender et al. 
199lh . 

/ t /max \ * / t /max \ 

V0-- e)/e (5) 


T/max 

*rot 


<Tn 


i/max 

'rot 
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where cr m is the mean velocity dispersion and e is the cllipticity of the cluster. Note that 
(Kot ax / cr m)* •C 1 indicate strong anisotropy. Using the relation above, we obtained the 
anisotropy parameter (V r j^ ax /cr m )* = 0.707 for the Ca-w group and 1.62 for the Ca-s group, 
where we adopted e = 0.074, V™ x = 1.0 km s -1 and o m = 5.0 km s -1 for the Ca-w group 
and e = 0.106, U r ” 1 t ax = 4.3 km s^ 1 and cr m = 7.7 km s -1 for the Ca-s group. For ellipticity, 
we used the simple mean of those given in Table [6] within about 210 arcsec from the center 
(i.e. about the half-light radius of the cluster). The large values of anisotropy parameters 
for both populations suggest that the flattening observed in M22 is most likely originated 
from the presence of internal rotation. 


5. SUMMARY 

We presented new wide-held ground-based Ca by photometry and the HST WFC3 
photometry of the peculiar GC M22 (NGC 6656), confirming our previous results that M22 
has a distinctive RGB split in the hk index, mainly due to the bimodal distribution in heavy 
elemental abundances, especially in calcium. 

Our ground-based photometry is complete down to V ~ 18.5 mag both in the central 
and the outer part of the cluster and incompleteness at the fainter magnitude regime does 
not affect our results presented in this paper. Although M22 is located toward the Galactic 
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bulge, the contamination from the off-cluster field populations appears to be negligible when 
our cluster membership star selection scheme based on our multi-color photometry is applied. 


Similar to our previous work (Lee et al. 2009_a), we defined the Ca-w and the Ca-s RGB 
groups based on the hk index values at a given V magnitude. We obtained the apparent 
number ratio of n(Ca-w):n(Ca-s) = 69:31 (± 3), without the correction for the differential 
evolutionary effects due to the heterogeneous elemental abundance s and slightly different 
ages between the two groups ( Marino et al. 2011 : Joo fe Lee 2013h . After applying such 
correction, the intrinsic or the true number ratio becomes Ji(Ca-w):n(Ca-s ) = 70:30 (± 3). 


Using the HST/WFC3 data of the central part of the cluster, we obtained the appar¬ 
ent SGB number ratio of n(bSGB):n(fSGB) = 69:31 (± 6) without the correction for the 
differential evolutionary effects and the true SGB number ratio of 62:38 (± 6) with the cor¬ 
rection, which is in excellen t agre ement with previous result from the HST observations in 
other passbands by Marino et ah (120091 ). 


We showed that the RGB split in the ml index in M22 is not mainly due to the difference 
in the CNO abundances, but in the heavy metal abundances. Remarkable resemblance in the 
HB morphology and the double RGB sequences both in the ml and the hk indexes between 
M22 and the composite GC using M55 and NGC 6752 (each of which has similar metal 
abundances as the Ca-w and the Ca-s groups in M22) may support our idea. Therefore, it 
is thought that the origin of the RGB split in the ml index in M22 is the same as that in 
the hk index: the difference in the mean heavy metal abundances. 


Comparisons with the chemical abundance measurements of RGB stars from previous 
studies showed that the hk index indeed a measure of calcium abundance at a given lumi¬ 
nosity. Our hk measurements are well correlated with Ca II H and K lines, infrared Ca II 
triplet lines and visual Ca I lines by others. Therefore, the hk index can provide an efficient 
means to distinguish the difference in the calcium abundance and furthermore the multiple 
populations in GCs and nearby dwarf spheroidal galaxies. 


Using homogeneous CN and CH band strength measurements by Smolinski et ah (20111). 
we showed that the gradients in <5CN and 5CH against metallicity are naturally expected to 
occur, when the common lower envelope for GCs with heterogeneous metallicity i s adopted. 


In this regard, we showed that the CN-CH positive correlation in M22 reported by Lim et al. 


( 201511 is most likely due to metallicity effects between the two RGB groups with heteroge¬ 
neous metallicities in M22. However, when corrected with its own lower envelope, each RGB 
group keeps maintaining its own CN-CH anti-correlation, showing a well-known character¬ 
istics of normal GCs. This may suggest that the merger of two GCs with different chemical 
compositions would be the most simplest explanation. 
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Our cy index measurements for each RGB group in M22 confirmed that there exists a 
significant spread in the nitrogen abundance among RGB stars in both groups. Quantita¬ 
tively, the cy index appears to over-estimate the nit roge n abunda nce w hen compared with 
the results from the high-resolution spectroscopy by Marino et al. (2011). 


We found that the mean RGB bump V magnitude, 14 U m P , of the Ca-w group is about 
0.15 mag brighter than that of the Ca-s group, which is consistent with the idea that the 
difference in Vb ump between the two groups is mainly due to the difference in metallicity. 
Furthermore, the tilted 14 ump magnitude level against A hk for the Ca-w group can be seen 
our data, suggesting that the hk index is indeed a good measure of metallicity and the RGB 
bump luminosity decreases with increasing metallicity. ft is an widely ac cepted f act t hat 


the formation of the EBffB stars may require a helium enhancement (e.g. D'A ntona et al. 


20021) . If the Ca-s group has an enhanced helium abundance by AY = 0.09 as Joo & Lee 


020131) proposed, the brighter 14 U m P magnitude by AV = 0.08 mag is expected for the Ca-s 
group. But this is not the case and the difference in the 14 ump magnitude within the context 
of the enhanced helium abundance of the Ca-s group is difficult to explain. 


With the carefully selected sample of RGB stars, we derived the coordinates of the center 
of both RGB populations using three different methods, finding insignificant differences in 
the coordinates of the center of both populations. The coordinates of the center of each 
group appear to be slightly different, but the difference is not substantially large to claim 
that the centers of two groups are distinctively different. On the other hand, the radial 
and spatial distributions between the two groups are different. The Ca-w group appears to 
be slightly more centrally concentrated while the number ratio between the Ca-w and the 
Ca-s groups are almost flat within the measurement errors up to ~ 3 — 4 half-mass radii, 
a reminiscence of the peculiar GC NGC 1851. The spatial distribution of the Ca-s group 
is very similar as that of the Ca-w group up to ~ 1.5 core radii, while the more elongated 
spatial distribution of the Ca-s group at larger radii can be seen. 


The similarity and the difference between the RGB and the ffB populations could be 
the critical puzzle to understand the formation of the EBffB stars. The BffB stars in M22 
are thought to be the canonical progeny of the metal-poor Ca-w RGB population with a 
normal helium abundance while the EBffB stars are though t to be tha t of the metal-rich 
Ca-s RGB population with an enhanced helium abundance ( Joo fe Leel 120131) . Our result 
showed that the BHB population is more centrally concentrated than the EBHB population, 
qualitatively consistent with the result from the RGB populations. However, the number 
ratio between the BHB and the EBHB populations is significantly different from that of the 
RGB and the SGB populations, n(BHB):n(EBHB) = 80:20 (± 10) without the correction 
for the differential evolutionary effects and n(BHB):n(EBHB) = 88:12 with such correction. 
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With this HB number ratio, we showed that the comparison of the modified R values of the 
two groups is against the theoretical expectation that the Ca-s + EBHB population (G2 
desi gnated by Joo fe Lee) is more helium enhanced than the Ca-w + BHB population (Gl) 


as 


Joo &; Lee ( 2Q13h proposed. 


Our results for the HB number ratio and the modified R values may indicate that 
the EBHB is not the only progeny of the C a-s RGB group. In this regard, the recent 


spectroscopic study of the HB stars in M22 by G ra t ton et al. (2 0141 ) is very intriguing as we 
have mentioned previously. They carried out the spectroscopic survey of the BHB stars in 
M22 and found a significant fraction ( ga 18%) o f the metal-rich helium-enhanced BHB stars 
in their sample (G r oup 3 designated bv lGratton et al.l) . If we consider the Group 3 BHB stars 
by Gratton et al. (2014) as the progeny of the Ca-s RGB stars, the apparent HB number 
ratio equivalent to the Ca-w and the Ca-s RGB groups becomes 66:34 and this number ratio 
between the two different HB groups is in good agreement with those from the RGB and the 
SGB starsj^l However, this does not explain why some of the metal-rich helium-enhanced 
RGB stars end up with the BHB stars, while others end up with the EBHB stars. Perhaps, 
may the discrepancy in the number counts indicate that there exists a non-canonical channel 
of th e evo lution of the_EBHB stars without invoking the helium enhancement as proposed 
by Moehler et al. (2 0111) ? Or, may this indicate that there exist the substantial spread or 
the bimodal distribution in helium abundance of the second generation of stars? Or, perhaps 
may the dispersion in the mass loss during the RGB phase be much larger than that of the 
canonical channel? These are the questions should be answered in the future. 


We measured the mean radial velocities and velocity dispersions of of each RGB pop¬ 
ulation using the results from a wide-held spectroscopic survey by ILane et alj (120091) . We 
obtained v r = —144.59 ± 0.40 km s -1 and a v = 4.71 ± 0.35 km s -1 for the Ca-w RGB 
population and v r = —145.77 ± 0.89 km s _1 and a v = 7.55 ± 0.66 km s” 1 for the Ca-s 
RGB population. The mean velocities of both populations are in good agreement, while the 
velocity dispersion of the Ca-s population is slightly larger. 


Finally, we derived the projected rotation velocities for each population using three 
different methods. Our rotation velocity measurement for all RGB stars is in good agreement 
with previous measurements by others, ranging from 1.4 km s -1 to 2.1 km s -1 depending 
on the methods, suggesting that the elongated structure of M22 is most likely due to the 
presence of internal rotation. However, detailed comparisons between the Ca-w and the Ca- 
s groups revealed that they have different kinematic properties. The Ca-s RGB population 


9 Recall that the Ca-s group in M22 may be equivalent to NGC 6752, which contains both the BHB and 
the EBHB stars. 
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(2.5 - 4.3 km s' 1 ) appears to rotate faster than the Ca-w population (1.0 - 1.5 km s' 1 ), 
with different position angles of the axis of the rotation, 89° for the Ca-w group and 64° for 
the Ca-s group. 


6. DISCUSSIONS 


The up-to-date observational evidences strongly suggest that the formation history of 
M22 must have been different from that of normal GCs in our Galaxy. Two different for¬ 
mation scenarios have been frequently explored in order to elucidate the nature of M22: the 


self-enrichment and the merger scenarios (e.g. iLee et al.l l2009al: iMarino et al.l 120 1 lL 12012 


Joo fe Lee 2013; Lim et al. 20151). Each hypothesis has pros and cons to explain the detailed 
photometric and spectroscopic characteristics observed in M22. 

In an attempt to delineate the chemical evolution between the two groups, at least 
three or four sources for specific elemental abundance changes should be involved: (i) Core- 
collapsed SNe: Heavy metals including Ca and Fe, but not affect the r-process elements, 
(ii) Intermediate-mass AGBs (IMAGBs) or FRMSs: Na-O, C-N anti-correlations. The He 
and overall C and N abundance enha ncements, (iii ) Low - mass. AG Bs (LMAGB): s-proces s 
elements. As extensively discussed bv IMarino et al.l 1)201 ll. 120121 1 and lRoederer et al.l (120111) . 
the self-enrichment formation hypothesis may require unlikely complex fine-tuning assump¬ 
tions. The formation scenario of M22 in the cont ext o f the self-enrichment hypothesis can 
be described in the following step0 (e.g Marino et ah 2011 1: (1) The formation of the first 
generation (FG) of stars of the s-process-poor (or Ca-w ) group (FGSP); (2) SNe II explo¬ 
sions, which expelled material far from the center; (3) IMAGBs of the FGSP polluted the 
intra-cluster medium, which went into the central region via cooling flow; (4) The formation 
of the second generation (SG) of stars of the s-poor group (SGSP); (5) LMAGBs of the FG 
expelled the ejecta with enhanced s-process elements; (6) A second cooling flow, containing 
the ejecta from the LMAGBs and SNell, into the central part of the cluster; (7) The for¬ 
mation of the FG of the s-process-rich (or Ca-s ) group (FGSR); (8) IMAGBs of the FGSR 
polluted the intra-cluster medium, which went into the central region via cooling flow again; 
(9) The formation of the SG of the s-rich group (SGSR); (10) Loss of stars only in the FGSP 
and the FGSR due to gas expulsion induced by SNe explosions. 

The first problem with such self-enrichment hypothesis is the number ratios between 
the FG and the SG in both the s-poor and s-rich groups. The most conspicuous feature 


10 The self-enrichment hypothesis by .Too fo Lee (201 3i) is slightly different from that of lMarino et ali J201ll ), 


but it does not affect our discussions made in this paper. 
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can be seen in the recent spectroscopic study of the majority of normal GCs is that t he SG 
is the major compo nent of normal GCs, generally ~ 50-80 % (e.g. iD’Antona et al.l 120081: 
Carretta et al.l 120091) . This cannot be easily explained with the chemical evolution with 


variant IMFs, which led D’ercole et ah 020 081 ) to propose the loss of the vast amount of the 
FG stars due to the gas expulsion induced by SNe la explosions in the early phase of the 
GC evolution. This must be the case for both groups of stars in M22. We estimate that the 
fractions of the SG for both gro ups in M22 ar e very similar, about 40-50%, based on the 
[O/Fe] versus [ Na/Fe) diagram of Ma rino et ah (201l|), using the definition for the primordial 
component by Carretta et ah ( 2009h . This means that the FGs in both groups are loosely 
bounded to the same extent with respect to each SG, so that the similar fractions of the FG 
stars in both groups must have been lost during the phase (10) above, which appears to be 
difficult to explain!^! 

The second problem is the number ratio between the two groups. Given the current 
theoretical framework of the normal GC formation which requires extensive mass loss of the 
previous generations of stars, it may be expected that the fraction of the s-rich (or the Ca-s ) 
stars is supposed to be larger than that of the s-poor (Ca-s) stars, in sharp contrast to the 
result presented here, n(Ca-w):n(Ca-s ) ~ 70:30. Initially, if M22 were much more massive 
and located in an isolated place, such as NGC 2419, it c ould b e ab le to retain the s-poor 
(or Ca-w ) stars expelled during the phase (10) above (e.g. Caloi & D’Antonal 20111). In fact, 
the recent o bservations showed that the metal-p oor population is the major component of 
NGC 2419 ( Beccari et al. 2013 : Lee et al. 2013 1. If so, isn’t it natural to expect that the 
FGSP and the FGSR are the major components for both groups? 

The third problem is the almost flat number ratio against the radial distance, with a 
hint of the slightly more centrally concentrated structure of the Ca-w group in M22. One of 
the general ch aracter istics of normal GCs is the central concentration of the SG population 

(e.g. 


Lardo et al. 2011 1. including NGC 2419 jBeccari et al. 2013 1. As mentioned earlier, the 


Wbody numerical simulations by Vesperini et ah ( 20131 ) suggested that the time required 
to have a flat number ratio between the FG and the SG of normal GCs is likely much larger 
than a Hubble time. If their result is applicable to the two groups of stars in M22, the flat 
number ratio against the radial distance is not likely due to the dynamical evolution_of the 


cluster with multiple stellar populations but it must have been set initially (e.g. Larsen et al. 
2015l h 


The last problem is the difference in the velocity dispersions. As discussed by Bejkki 


n The loss of the FGSP can occur during the phases (4-5), but this cannot explain the similar fractions of 
the FG in both groups and the number ratio between the two groups. 
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(120101 ). the SG formed gaseous ejecta of the FG in the central region of the cluster can have 
larger rotational amplitudes but smaller velocity dispersions than the FG. He suggested that 
the dynamically cold nature of the SG is largely due to gaseous dissipation during the gas 
accretion onto the cluster center. Thus the lower velocity dispersion of the SG is the key 
signature that can be imprinted. If the Ca-s group were formed out of gaseous material 
accreted onto the center within the self-enrichment scenario, the lower velocity dispersion is 
expected. However this is not consistent with our result, showing that the velocity dispersion 
of the Ca-s group is slightly larger than that of the Ca-w group. Therefore, the difference in 
the velocity dispersions between the two groups does not appear to be well understood within 
the current frame of the dynamical evolution of the GCs with multiple stellar populations. 

It is interesting to note that the heavy elemental abundance pattern between the two 
groups of M22 can be reproduced by the merger of two mono-metallic normal GCs, M5 
and M4 (IRoederer et all 2011b In addition, we showed that the photometric characteristics 
of M22 can also be successfully reproduced by the two mono-metallic normal GCs, M55 
and NGC 6752. As mentioned earlier, the observed CN and CH distributions in M22 RGB 
stars can be simply explained with the merger scenario. The differences in the spatial 
distributions and the kinematic properties between the two groups in M22 can not be well 
understood in the context of the self-enrichment hypothesis and they appear to favor the 
merger scenario. Also, the rather high flattened structure and rotation of M22 can be induced 
by the merging of a GC pair. We proposed that the plausible explanation of the formation 
of M22 would be that M22 formed via a merger of two GCs, preferentially.in a dwarf galaxy 


environment where the G C merger events are much more probable (e.g. iThurl fc Johnston 


2002 


2010 


Marino et ah 


20U). 


Bekki fe Yondl2012h . and then accreted later to the Milky Way (see also lCarretta et al 
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Table 1. Journal of observations. 


date 

(yy/mm/dd) 

Int. Exp. Time (s) 
y b Ca 

Tel. 

(CTIO) 

Weather 

11/03/27 

180 

310 

1,500 

1.0m 

Photometric 

11/03/28 

340 

580 

2,400 

1.0m 

Non-photometric 

11/03/29 

340 

580 

2,400 

1.0m 

Photometric 

11/03/30 

340 

580 

1,800 

1.0m 

Non-photometric 

11/03/31 

160 

540 

2,400 

1.0m 

Photometric 

11/04/01 

600 

1,350 

4,200 

1.0m 

Non-photometric 

11/04/02 

275 

600 

2,100 

1.0m 

Photometric 

11/04/03 

400 

900 

3,000 

1.0m 

Photometric 

11/04/07 

200 

450 

1,500 

1.0m 

Photometric 

11/08/19 

100 

200 

1,200 

1.0m 

Photometric 

11/08/20 

200 

400 

1,500 

1.0m 

Non-photometric 

11/08/21 

350 

600 

2,400 

1.0m 

Photometric 

11/08/22 

300 

500 

2,400 

1.0m 

Photometric 

11/08/31 

310 

240 

900 

1.0m 

Non-photometric 

11/09/04 

100 

200 

900 

1.0m 

Photometric 

12/04/29 

60 

180 

1,200 

0.9m 

Non-photometric 

12/07/12 

250 

330 

1,200 

0.9m 

Non-photometric 

12/07/14 

200 

600 

1,200 

0.9m 

Non-photometric 

12/07/16 

100 

300 

1,500 

0.9m 

Non-photometric 

12/07/18 

200 

500 

1,500 

0.9m 

Non-photometric 

12/07/19 

100 

250 

1,500 

0.9m 

Non-photometric 

12/07/21 

200 

500 

1,500 

0.9m 

Non-photometric 

12/07/22 

60 

180 

2,400 

0.9m 

Non-photometric 

12/07/23 

50 

150 

1,200 

0.9m 

Photometric 

12/07/24 

130 

390 

1,200 

0.9m 

Photometric 

12/10/10 

140 

420 

3,000 

1.0m 

Photometric 

12/10/11 

120 

360 

2,400 

1.0m 

Photometric 

13/04/07 

80 

200 

1,200 

1.0m 

Photometric 

13/04/08 

185 

475 

1,500 

1.0m 

Non-photometric 

13/04/11 

180 

500 

1,500 

1.0m 

Photometric 
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Table 1—Continued 


date 

(yy/mm/dd) 

Int. Exp. Time (s) 
y b Ca 

Tel. 

(CTIO) 

Weather 

13/04/12 

115 

310 

1,500 

1.0m 

Photometric 

13/04/13 

100 

250 

1,200 

1.0m 

Photometric 

13/05/14 

100 

250 

1,500 

1.0m 

Non-photometric 

13/07/27 

80 

200 

1,200 

1.0m 

Non-photometric 

13/07/28 

180 

260 

1,200 

1.0m 

Non-photometric 

13/07/31 

95 

240 

1,200 

1.0m 

Non-photometric 

13/08/02 

180 

500 

1,500 

1.0m 

Non-photometric 

14/05/23 

150 

450 

1,800 

1.0m 

Non-photometric 

14/05/25 

300 

900 

3,900 

1.0m 

Non-photometric 

14/05/26 

300 

900 

3,900 

1.0m 

Photometric 

14/05/27 

240 

600 

3,600 

1.0m 

Photometric 

14/05/30 

320 

800 

3,600 

1.0m 

Photometric 

14/05/31 

615 

1,545 

8,700 

1.0m 

Photometric 

14/06/01 

120 

360 

1,500 

1.0m 

Photometric 

total 

9,145 

20,930 

90,900 
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Table 2. Comparisons of photometric results. 



Anthony-Twarog et al. 

Richter et al. 

No. of stars 

198 

60 

AV 

+0.024 ±0.051 

-0.038 ±0.027 

A(b-y) 

±0.034 ±0.050 

±0.003 ±0.037 

Ami 

-0.028 ±0.063 

±0.021 ±0.037 

A hk 

-0.029 ±0.079 
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Table 3. Color-Magnitude Diagram data. 


ID 

RA 

Dec 

V 

b-y 

hk 

96 

279.502014 

-24.007305 

11.436 

0.849 

1.280 

104 

279.061249 

-24.014750 

11.895 

0.717 

0.901 

106 

279.059052 

-24.033777 

11.756 

0.927 

0.954 

130 

278.959198 

-24.261583 

11.649 

0.895 

1.178 

146 

279.027954 

-24.347834 

12.278 

0.528 

0.463 

163 

279.502747 

-24.352249 

11.570 

0.836 

1.186 

197 

279.064087 

-23.805361 

12.211 

0.885 

0.815 

252 

279.312927 

-23.539223 

11.941 

0.308 

0.409 

285 

279.184082 

-23.980139 

12.530 

0.888 

0.842 

297 

279.138306 

-23.940722 

12.655 

0.942 

0.830 

301 

278.767181 

-24.275888 

13.256 

0.202 

0.221 


Note. — Table [3] is presented in its entirety in the 
electronic edition of the Astrophysical Journal Supple¬ 
ment. A portion is shown here for guidance regarding 
its form and content. 
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Tabic 4. Number ratios. 



Before Correction 

After Correction 

Error 

n(Ca-w):n(Ca-s) 

69:31 

70:30 

± 3 

n(bSGB) :n(fSGB) 

69:31 

62:38 

± 6 

n(BHB) :n(EBHB) 

80:20 

88:12 

± 10 
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Table 5. Centers of M22. 



Simple mean 

Half-sphere 

Pie-slice 


AX(") 

AY (") 

AX(") 

AY(") 

AX(") 

AY (") 

All 

2.9 

-0.1 

5.5 

1.1 

7.1 

1.0 

Ca-w 

3.8 

0.4 

7.5 

1.6 

7.7 

-2.8 

Ca-s 

0.8 

5.6 

-1.4 

-1.4 

-4.9 

6.5 
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Table 6. Ellipse fitting parameters. 


Pop. 

grid 

Ellipse center 

9 

b/a 

e 



AX(") 

AY(") 




Ca-w 

0.9 

10.5 ± 0.4 

4.2 ± 0.4 

1.6 ± 1.4 

0.929 ± 0.016 

0.071 


0.8 

10.1 ± 0.6 

4.6 ± 0.6 

3.2 ± 1.3 

0.931 ± 0.015 

0.069 


0.7 

9.2 ± 0.7 

5.2 ± 0.7 

4.9 ± 1.2 

0.937 ± 0.014 

0.063 


0.6 

7.9 ± 0.7 

5.0 ± 0.8 

7.3 ± 1.1 

0.940 ± 0.013 

0.060 


0.5 

6.2 ± 0.8 

4.3 ± 0.8 

11.2 ± 1.0 

0.943 ± 0.011 

0.057 


0.4 

3.5 ± 0.8 

2.7 ± 0.8 

17.4 ± 0.8 

0.942 ± 0.010 

0.058 


0.3 

0.2 ± 0.7 

-0.6 ± 0.8 

24.3 ± 0.7 

0.934 ± 0.008 

0.066 


0.2 

-2.3 ± 0.7 

-5.2 ± 0.7 

31.4 ± 0.5 

0.912 ± 0.006 

0.088 

Ca-s 

0.9 

-6.6 ± 0.4 

6.2 ± 0.4 

15.3 ± 1.3 

0.943 ± 0.015 

0.057 


0.8 

-3.3 ± 0.6 

5.4 ± 0.6 

18.9 ± 1.2 

0.943 ± 0.014 

0.057 


0.7 

0.1 ± 0.7 

3.6 ± 0.7 

24.1 ± 1.1 

0.944 ± 0.013 

0.056 


0.6 

3.8 ± 0.8 

0.8 ± 0.7 

29.2 ± 1.0 

0.941 ± 0.011 

0.059 


0.5 

6.8 ± 0.8 

-3.2 ± 0.8 

35.4 ± 0.8 

0.934 ± 0.010 

0.066 


0.4 

8.1 ± 0.8 

-8.0 ± 0.8 

42.7 ± 0.7 

0.912 ± 0.008 

0.088 


0.3 

6.6 ± 0.7 

-13.2 ± 0.8 

41.8 ± 0.6 

0.872 ± 0.006 

0.128 


0.2 

2.4 ± 0.7 

-13.5 ± 0.7 

37.6 ± 0.4 

0.837 ± 0.005 

0.163 
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Fig. 


fa) A co mparis on of Ca filter transmission functions between that 


m 


Anthony-Twaro g et al. (119911) (the dotted line) and that for our new Ca filter measured 


with collimated beam (the solid line). Both filters have similar FWHMs, approximately 90 
A but our new Ca filter has a more uniform and high transmission across the passband, 
dropping more rapidly at both edges. The transmission function for the new filter system, 
JWL39, is also shown in the plot, (b) Synthetic spectra for the CN normal (the blue line) 
and the CN strong (the red line) RGB stars. The CN band at A 3885 A lies outside of the 
lower boundary of our new Ca filter. On the other hand, we can measure the CN band 
strength by using our new filter system, JWL39. (c) Filter transmission functions for the 
ground-based extended Stromgren photometry, (d) Filter transmission functions for the 
HST WFC3/UIVS photometry. Note that the bandwidth of the HST F547M is about twice 
as broad as that of the Stromgren y filter. 
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Fig. 2.— Residu als of the CCD photometry of lAnthonv-Twarog et al.l (119951) and 


Ri chte r et al . (119991) as a function of the magnitude. The differences are in the sense of 


Anthony-Twarog et al. or Richter et al. minus our work. 
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Fig. 3.— Completeness fractions at each magnitude bin measured from artificial star exper¬ 
iments for the inner and outer parts of M22. The experiments suggest that our photometry 
is complete down to V ~ 18.5 mag. 
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Fig. 4.— Color-magnitude diagrams for bright stars located within 10 arcmin from the 
center of the cluster. Note that discrete two RGB sequences can be clearl y seen in the hk 
versus V CMD. In the figure, the ml versus V CMD is taken from Lee et al. (2009a ). 
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Fig. 5.— (a) & (b) CMDs of the M22 field. The enclosed regions by blue lines indicate the 
location of RGB stars in M22, while the enclosed region by cyan box in (a) mostly represents 
the bulge red-clump star population; (c) The number of stars in the enclosed region is 3571. 
Compared with Lane e t al. 020091 1. we have 307 stars in common; 295 stars are the radial 
velocity member star (green dots) and 12 are radial velocity non-member star (red plus 
signs), (d) The hk versus V CMD for 3571 stars from (c). The number of stars in the 
enclosed region is 2054 without the radial-velocity non-member stars. 
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1 2 3 

log (r/arcsec) 


Fig. 6.— The surface-brightness profiles for M22 RGB stars with 12.0 < V < 16.5 mag. 
(a) Ca-w RGB stars only; (b) Ca-s RGB stars only; and (c) All RGB stars in M22. Our 
surface-brightness profiles for M22 RGB stars are in excellent agreement with the Chebyshev 
polynomial fit of M22 bv iTrager et al.l (119951 1. shown with green solid lines. On the other 
hand, the surface-brightness profile for the Galactic bulge red-clump stars is almost flat 
against the radial distance from the center of the cluster. Our results strongly suggest that 
the contamination from the off-cluster field star in our RGB selection procedure is not severe. 
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Ahk 


Fig. 7.— (Left) The plot of Ahk versus V for M22 RGB stars, where Ahk is defined to 
be the difference in color from the mean fiducial sequence of the Ca-w group. Stars with 
the probability of being the Ca-w group P(w) > 0.5 from the EM estimator assuming the 
two gaussian mixture model are denoted with blue dots, while P(s) > 0.5 with red dots. 
(Right) Histograms of M22 RGB stars with different magnitude bins. The two gaussian 
distribution was estimated by using the EM method and the number ratio in the parenthesis 
is for n(Ca-w):n(Ca-s) in each bin. 

































- 2.0 - 1.0 0.0 1.0 2.0 3.0 4.0 


My 

Fig. 8.— The theoretical luminosity functions for the Ca-w group (the blue line) and the 
Ca-s group (the red line). The grey shaded area corresponds to the RGB stars and the yellow 
shaded area to the SGB stars analyzed in our current work. Note that the RGB number 
ratio between the two groups at the fixed magnitude bin is very similar, n(Ca-w):n(Ca-s) 
= 48.3:51.7 (± 0.3), while the number of SGB stars in the Ca-w (or bSGB) group is slightly 
larger than that in the Ca-s (or fSGB) group, n(bSGB):n(fSGB) = 58.0:42.0 (± 0.2). 
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Fig. 9.— (Upper panels) CMDs of M22 using stars with 0.5 < r < 5.0 arcmin. In the ml 
versus V CMD, we show stars with (b — y) > 0.45 mag, since the location of HB stars in 
the CMD overlaps with that of RGB stars. Note that the blue_solid_line is foi_the Padova 
isochrone with [Fe/H] = —1.8 and the red solid line for —1.5 (Bre ssa n et ah 2012). (Bottom 
panels) Composite CMDs for globular clusters M55, which mimics the Ca-w groups in M22, 
shown with blue dots and NGC 6752, which mimics the Ca-s groups, with red dots. The 
numbers in the parentheses indicate differential colors and the visual magnitude terms for 
NGC 6752 with respect to M55. The composite CMDs using M55 and NGC 6752 reproduce 
very well the photometric characteristics of M22 in the ml versus V and hk versus V CMDs. 
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Fig. 10.— Color-magnitude diagrams of M22 using the HST data. The /i/cstmag denotes 
(?uf 395 n — ?71f467m) — ('W-F 467 M — W-F 547 M) • The split in the SGB and the RGB sequences in 
M22 can be clearly seen in the /i/cstmag versus m-p 547 M CMD. 
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Fig. 11.— (Left panels) Artificial CMDs around the SGB region of M22 for different pho¬ 
tometric systems. The blue dots are for the bSGB stars and the red dots are for the fSGB 
stars. The vertical and horizontal grey shaded regions indicate ranges in the color indexes 
and in the magnitudes used to derive the number ratios between the two groups. We show 
only a small fraction of stars in each plot for the sake of clarity. (Middle panels) The dis¬ 
tributions of SGB stars at the fixed color indexes. The first number ratios in each plot 
are those returned from the EM estimator assuming the two gaussian mixture model while 
the number ratios in the parenthesis are the input values from our evolutionary population 
models. (Right panels) The distributions of SGB stars at the fixed magnitudes. Note that 
the number ratios from the ( hk, V ) plane are the most reliable, while those from the (b — y , 
V) plane are the least reliable. 
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Fig. 12.— (a) The hfcsTMAG versus nip 547 m CMD around the SGB region of M22. Also 
shown is the fiducial sequence for the bSGB group. We use the SGB stars with 17.0 mag < 
m F 547 M < 17.4 mag to estimate the number ratio between the two groups. The stars with 
P(bSGB) > 0.5 from the EM estimator are denoted with blue dots, while P(fSGB) > 0.5 
with red dots, (b) The plot of A/i/c STM ag versus V, where A/i/c STM ag is defined to be the 
difference in the hk STMAG index between the individual stars and the fiducial sequence of the 
bSGB group, (c) The distribution of A/i/cstmag and the two gaussian mixture model fit to 
the data. We obtain the apparent number ratio of n(bSGB):n(fSGB) = 69:31 (± 6) and the 
intrinsic number ratio of n(bSGB):n(fSGB) = 62:38 (± 6) after applying the correction for 
the differential evolutionary effect. This SGB number ratio is in agreement with the RGB 
number ratio within measurement errors. 
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Fig. 13.— A comparison of the Calcium index A(Ca) measured by Norris & Freeman ( 19831 ) 
between the RGB stars in the Ca-w (blue plus signs) and the Ca-s (red circles) groups in 
M22 as functions of V magnitude. We also show the least square fit to the Ca-w RGB 
stars and the differences are measured from the fitted line to correct the temperature and 
the surface gravity effects on the absorption strengths. The histograms for each group are 
shown in the right panel. Note that the boundary between the two groups is rather sharp. 
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Fi&. 14. — Comparis ons of Ca II triplet strength E(CaT), (= W 854r2 + IT 8662 ), measured 
by Da Costa et ah ( 20091 ) for the Ca-w (blue circles) and the Ca-s (red circles) RGB stars 
against magnitude difference from the horizontal branch, V — Vhb- In the upper panel, the 
dotted line denotes the least square fit to the data. The lower panel show the residuals 
around the fitted line. 
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Fig. 15.— Plots of [Ca/H] versus [Fe/H] and [Na/Fe] versus [O/Fe] for the Ca-w (blue pluses) 
and the Ca-s (r ed s quares) RGB stars using high-resolution spectroscopic measurements by 
Marino et al. ( 2011 1. The plot of [Ca/H] versus [Fe/H] shows that the RGB stars in the 
Ca-w group have lower calcium and iron abundances than those in the Ca-s group. The 
plot of [Na/Fe] versus [O/Fe] shows that each group has its own Na-0 anti-correlations. In 
the mean the Ca-w group appears to have a higher oxygen and a lower sodium abundances 
than the Ca-s group. 
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Fig. _16.— Distributions of CN(3839) and CH(4300) against V — Vhb for 8 GCs by 


(ISmolinski et al.l l201ll) . The numbers inside the pare ntheses are the m etallicity of each 
GC. Also shown are two GCs, M22 and NGC 288, by Tam et ah ( 2015 ). The thick solid 
lines denote the common lower envelopes for RGB stars around Vhb, — 1 < V — Vhb <0.5, 


used i n our analysis. Note that the CN(3839) and_CH(4300) index values bv ISmolinski et al. 
(2011) do not agree with those by Lim ct ah (12015 ) due to different definitions of the band 
strengths between the two works. 
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Fi g. 17 .— The dist r ibutio ns of the <5CN(3839) and the hCH(4300) against [Fe/H] for 8 GCs 
by Smolinski et ah ( 2011 ). where the hCN(3839) and the hCH(4300) are defined to be the 
differences in the band strengths from the common lower envelopes. The least square fits to 
the data are shown with dashed lines. As expected from Figure [I6l both the hCN(3839) and 
the <5CH(4300) have the gradients against metallicity. 
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Fig. 18-— Distributions of CN(3839) and CH(4300) against V — Vhb for M22 RGB stars 
by Lim ct ah (120151 ). The blue color and the red color denote the Ca-w and the Ca-s RGB 
stars in M22, respectively. The blue dotted lines denote the lower envelope for the Ca-w 
stars, while the red dotted lines for the Ca-s stars. The vertical dashed lines denote the 
region around the Vhb of M22, — 1 < V — Vhb < 0.5. 
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Fig. 19.— The distributions of the hCN(3839) versus the hCH(4300); (a) Using the common 
lower envelopes, (b) The Ca-w RGB stars using its own lower envelopes, (c) The Ca-s RGB 
stars using its own lower envelopes. The RGB stars around the Uhb are shown with the 
plus signs and crosses. The hCN-hCH positive correlation in (a) is most likely due to the 
metallicity effect between the two groups as shown in Figure fT71 
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Fig. 20.— Plots of cy versus V CMDs of M22 RGB stars. The V magnitude of the horizontal 
branch, Vhb = 14.15 mag, is indicated by an arrow. The cy spread at Vhb is cr(cy) ~ 0.03 
mag in each group, corresponding to er[N/Fe] « 0.5 dex, if the spread in the NH absorption 
strengths are solely responsible for the cy spreads in M22 RGB stars. In the right panel of 
the figure, we show the cumulative LFs for the Ca-w and the Ca-s populations. We show 
the V magnitude levels of the RGB bump, V^, umpi at which the slope of the LF changes 
abruptly, for each population by arrows. We obtained Vb ump = 13.91 mag for the Ca-w 
population while 14.06 mag for the Ca-s population. The fainter bump magnitude in the 
Ca-s population suggests the metallicity effect on the bump brightness is greater than the 
effect from the helium enrichment is. 
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Fig. 21.— Generalized histograms for RGB stars against V magnitude, (a) Ca-w RGB stars 
with A hk < 0 from Figured (b) Ca-w RGB stars with A hk > 0. (c) Ca-s RGB stars. The 
vertical thin lines represent V\, um p f° r eac h subset, 13.90, 13.96 and 14.06 mag respectively. 
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Fig. 22.— Relative positions of centers with respect to that of Goldsbur v et al.1 (120101 ). The 
blue color is for the Ca-w group, the red color for the Ca-s group and the black for the 
combined RGB stars. Diamond signs denote for the centers from the simple mean, plus 
signs for those from the half-sphere method and crosses for those from the pie-slice method. 
The filled circles are centers of each population from the iso-density contour method from 
90% to 30% levels. 
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Fig. 23.— (a)The number ratios of the Ca-s to the Ca-w groups against the radial distance 
from the center. The RGB number ratio appears to have a weak radial gradient up to r « 
700 arcsec from the center, in the sense that the number of the Ca-w population is slightly 
more centrally concentrated than that of the Ca-s population. In the figure, the vertical 
lines denote the cluster’s core and half-light radii and the horizontal lines denote the mean 
number ratio and the standard deviations, (b) Cumulative radial distributions of Ca-w 
(blue) and Ca-s (red) RGB populations. The Ca-w RGB stars are slightly more centrally 
concentrated. The K-S test indicates a probability of 0.1% that two populations are drawn 
from the same parent population. 
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Fig. 24.— Spatial distributions of the Ca-w and the Ca-s RGB stars in M22 are shown in 
the upper panels and the smoothed contour maps using the fixed Gaussian kernel are shown 
in the lower panels, where we show the iso-density contour lines for 90, 70, 50, and 30% of 
the peak values for both populations. We also show the FWHM of our adopted Gaussian 
kernel in the lower left panel of the figure. Note that the distribution of the Ca-s RGB stars 
is more elongated. 
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Fig. 25.— The run of the axial ratio, b/a, the ellipticity, e (= 1 — b/a) of the Ca-w (blue) 
and the Ca-s (red) groups against the semi-major axis, a. The radial distributions for each 
group bifurcate at the radial distance larger than a > 150 arcsec, in the sense that the Ca-s 
group is more elongated. 
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Fig. 26.— CMDs of the HB region of M22 and the cumulative LF of HB stars. The HB 
gap at V = 15.97 mag, with the zero slope in the cumulative LF of HB stars, is denoted 
by arrows. The BHB and EBHB stars are denoted by blue filled circles and red plus signs, 
respectively. 
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Fig. 27.— Same as Figure [23j but for the BHB and EBHB stars in M22. Similar to the 
Ca-w RGB population, the BHB population is more centrally concentrated. 
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Fig. 28.— (Upper) The distribution of radial velocities of individual stars against the radial 
distance from the center.The blue color represents the Ca-w group and the red color the 
Ca-s group. No radial gradient can be seen. (Lower) Histograms for the radial velocities of 
each group with gaussian fits. The mean radial velocities and the velocity dispersions are 
very similar between the two groups. 
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Fig. 29.— Radial velocities of individual RGB stars in M22 against their position angles. 
The two numbers in the parentheses in each plot denote the mean radial velocity and the 
amplitude of the mean rotational velocity. The red solid lines are sinusoidal fits to the data. 
fa) All RGB stars; (b) Ca-w stars; and (c) Ca-s stars with radial velocities measured by 
Peterson fe Cudworth (119941) . (e) - (f) As (a) - (c) but using stars measured bv lLane et al. 


(120091 ). It is evident that the Ca-s group rotates faster than the Ca-w group does. 
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Fig. 30.— The difference in the mean radial velocities between both hemispheres against 
its position angle. Again, it is evident that the Ca-s group appear to rotate faster than the 
Ca-w group does. The position angle refers to that of the equator of the rotation, i.e. the 
axis perpendicular to the rotation axis. The error bar is the error of the mean. 
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Fig. 31.— The upper panels show the distributions of RGB stars with radial velocity mea¬ 
sured by Lane et ah! (120091) . The blue color denotes the blue shift and the red color the 
red shift. The size of each circle indicates the difference between the velocities of individual 
stars and the mean velocity. The bottom panels show the radial velocities and the velocity 
dispersions against the projected distance on the equator of the rotation. In the Figure, the 
position angles of 86° and 73° for the Ca-w and the Ca-s groups, respectively, are used. In 
the middle panels, red dotted lines represent the mean rotation. The maximum rotation 
velocities at the half-light radius are 1.0 km s -1 for the Ca-w and 4.3 km s -1 for the Ca-s 
groups. The velocity dispersion of the Ca-s group is slightly larger but those of both groups 
are in agreement within the errors. 



















































